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In QCD with two flavours of massless quarks, the chiral phase transition is plan- 
ch I sibly in the same universahty class as the classical four component Heisenberg 

magnet. Therefore, renormalization group techniques developed in the study of 
phase transitions can be applied to calculate the critical exponents which char- 
ts ■ acterize the scaling behaviour of universal quantities near the critical point. As 

a consequence of this observation, a quantitative description of the universal 
physics of the chiral phase transition in circumstances where the plasma is close 
I to thermal equilibrium as it passes through the critical temperature can be ob- 

tained. This approach to the QCD phase transition has implications both for 
CO ' lattice gauge theory and for heavy ion collisions. Future lattice simulations with 

^ . longer correlation lengths will provide quantitative measurements of the various 

exponents and the equation of state for the order parameter as a function of 
. temperature and quark masses. Present lattice simulations already allow many 

qualitative tests. In a heavy ion collision, the consequence of a long correlation 
length would be large fluctuations in the number ratio of neutral to charged 
pions. Unfortunately, we will see that because of the explicit chiral symmetry 
breaking introduced by the masses of the up and down quarks, no equilibrium 
(— I ' correlation length gets long enough for this phenomenon to occur. In the third 

chapter of this article, I discuss attempts to model the dynamics of the chiral 
order parameter in a far from equilibrium QCD phase transition by consider- 
ing quenching in the 0(4) linear sigma model. I argue, and present numerical 
^ ' evidence, that in the period immediately following the quench long wavelength 

modes of the pion fleld are amplified. This could have dramatic phenomeno- 
logical consequences in heavy ion collisions. I review various recent theoretical 
advances — attempts to include expansion and relax the quench assumption; at- 
tempts to include quantum mechanical effects. Long wavelength pion oscillations 
have now been seen in a number of simulations in which different assumptions 
are made. However, all of the theoretical treatments involve idealizations and 
assumptions, and should only be regarded as qualitative guides to what can hap- 
pen. It is up to experimentalists to determine whether or not such phenomena 
occur; detection in a heavy ion collision would imply an out of equilibrium chiral 
transition. 
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1. Introduction 

Shortly after the discovery that QCD is asymptotically freeB and that therefore 
quarks are almost noninteracting at short distances, Collins and Perryi noted that this 
means that QCD at high temperature or high density is qualitatively different than 
at low temperature and density. At high temperatures (T ^ Aqcd) or high baryon 
number density (n ^qcd) QCD describes a world of weakly interacting quarks 
and gluons very different from the hadronic world in which we live. This raises the 
possibility of phase transitions as the temperature or density is increased. While 
the subject of phase transitions as a function of increasing density has interesting 
applications in the physics of neutron stars and heavy ion collisions, in this article I 
will limit myself to discussing the physics of a QCD plasma with zero baryon number 
density. There are at least two qualitative differences between such a plasma at 
T 3> Aqcd and at T ~ 0. First, as Collins and Perry discussed, at low temperatures 
one has a plasma of quarks and gluons confined in mesons and baryons, while at high 
temperatures the quarks and gluons are deconfined. Second, at low temperatures a 
qq condensate spontaneously breaks chiral symmetry while at high temperatures the 
interactions among the quarks and anti-quarks are weak, no such condensate exists, 
and chiral symmetry is manifest. Thus, in studying the phase structure of QCD 
as a function of increasing temperature, there are two possible phase transitions to 
consider. 

Studying the QCD phase transition(s) has a certain intrinsic appeal. The question 
involved is simple enough that it can be stated in a sentence: What happens if you 
take a chunk of any substance around you and heat it up until it is one hundred 
thousand times hotter than the centre of the sun? While the question is simple, 
if one tries to answer it experimentally by colliding heavy ions at ultra-relativistic 
energies, the phenomena seem at first glance complicated and intractable. And yet, 
using techniques from different areas of theoretical physics — quantum field theory, 
particle physics, nuclear physics, lattice gauge theory, condensed matter physics, and 
cosmology — much progress is possible. Quantitative predictions for processes in 
which a plasma cools through the QCD phase transition in thermal equilibrium can 
be extracted, and striking phenomena with qualitative signatures may occur if the 
plasma is far from equilibrium as it cools through the transition. 

In addition to having intrinsic appeal, the QCD phase transition is of interest from 
several different points of view. First, there can be no doubt that it occurred through- 
out the universe about a microsecond after the big bang. Second, it is reasonable to 
hope that in a heavy ion collision of sufficiently high energy, a small region of the high 
temperature phase is created which then cools through the phase transition. Third, 
lattice gauge theory is well suited to calculating the equilibrium properties of QCD 
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at high temperatures. From all these perspectives, it is important to learn as much 
as can be learned analytically about the phase transition, relying as much as possible 
only on fundamental symmetries and universality arguments and as little as possible 
on specific assumptions and models. In Chapter 2, following original work done with 
Wilczek,i we will see that it is possible to use this approach to learn much about the 
equilibrium QCD phase transition. In studying heavy ion collisions, however, it is 
worth considering the possibility that the transition process may be far from thermal 
equilibrium. Without thermal equilibrium, it is much harder to be quantitative both 
because the phenomena are more varied and because there is no notion of universality. 
We will nevertheless find that the distinctive qualitative signature of the QCD phase 
transition in heavy ion collisions discussed in Chapter 3 — the amplification of long 
wavelength oscillations of the pion field — can only occur if the plasma is far from 
thermal equilibrium as it experiences the QCD phase transition.! 

The physics of QCD at high temperatures is an enormous subject with many 
facets, which I will not attempt to review here. The reader interested in entering 
the voluminous literature might consider starting with the review by Gross, Pisarski, 
and Yaffe,i the review by McLerranJ§ the book by Shuryak,lll the previous collection 
of articles edited by Hwa,i recent reviews by Miiller,! recent conference proceedings 
in the Quark Matter series,0inilii or other articles in this volume. 

The remainder of this introduction has two goals. First, a brief discussion of 
the confinement /deconfinement phase transition seems in order, since the bulk of 
this article is concerned with the chiral phase transition. Second, I will give a very 
brief introduction to some of the physics relevant for a discussion of the QCD phase 
transition in heavy ion coUisons, in order to set the stage for the third chapter. 

Let us recall the arguments^ which lead to the conclusion that the confine- 
ment /deconfinement phase transition is first order in pure SU{3) gauge theory and is 
smoothed into a crossover when light quarks are present. In studying SU{N) gauge 
theory in thermal equilibrium at a temperature T, the partition function is given by 
a Euclidian path integral 

Z = [ VA^exp{~SE) (1) 

JPBC 

over a slab of spacetime which is infinite in the three spatial directions and is /3 = 1/T 
thick in the t direction. = AJJA*^ where the A° are representation matrices for 
the fundamental representation of SU{N). The path integral is over gauge field 
configurations which satisfy periodic boundary conditions A^{x,0) = A^{x,(3). The 
Euclidian action is given by 

SE = ^[^tj d'f trF^^F^, . (2) 
The expectation value of operators is given by 

{0{A^)) = Z-^[ VA^eM-SE)0{A^) . (3) 

JPBC 
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We will be particularly interested in the Wilson loop operator 

r ff^ 1 

L{x) = tr Q{x) ; Q{x) = Pexp i / dtAo{x,t) (4) 

L JO J 

for a loop at fixed x which loops around the periodic time direction once. P denotes 
path ordering. This operator, conventionally called the Polyakov loop, is of interest 
because the free energy of a single static quark is 

/3(F,-Fo) = -log(L(x)) (5) 

where Fq = —TlogZ is the free energy of the vacuum. If (L) = 0, inserting a quark 
into the ensemble costs an infinite free energy. This implies that the flux cannot be 
screened — the state of the fields is changed all the way out to spatial infinity 
and costs an infinite free energy. Therefore, (L) = indicates that the theory is in 
the confining phase. If (L) ^ 0, then Fg is finite and the theory is in a deconfined 
phase characterized by Debye screening. 
Under a gauge transformation 

UA^U-' + lUd^U-^ , f/(f , t) e SU{N) (6) 

L transforms as 

L{x) trf/(f , 0)Q{x)U{x, p) , (7) 

so L is invariant when U is periodic. Local observables are invariant under those 
aperiodic gauge transformations 

U{x,t + p) = zU{x,t) (8) 

where z = exp{2TTin/N) is an element of the Z{N) centre of the gauge group. This 
is not true of the Polyakov loop, which transforms under according to 

L^zL . (9) 

Therefore, (L) is an order parameter. Going from (L) = to (L) ^ requires the 
spontaneous breaking of a symmetry — in this case, spontaneous breaking of Z{N). 
Hence, the realization of the Z{N) symmetry directly determines whether or not the 
SU (N) gauge theory is in a confining phase. 

The confinement /deconfinement phase transition as a function of temperature in 
SU {N) gauge theory is a transition in which a global Z{N) symmetry is spontaneously 
broken. If the order parameter (L) varies continuously (but non-analytically) from 
(L) = for T < Tc to (L) ^ for T > Tc, then the transition is second order. If, 
on the other hand, (L) is discontinuous at T^^ the transition is first order. Following 
Landaujli we write an effective Lagrangian for L at long length scales whose terms 
are determined by the symmetry of the order parameter 

^eff = / d^x\{^^^f + (10) 
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where the effective potential is 

ViL) = fi\L\')+gi\L\^Re{L'^)) . (11) 

In general, / and g are arbitrary functions. To look for a second order phase transition, 
we expand them about the origin giving 

V{L) = a|L|2 + b\L\^ + cRe{L^) + ■■■ (12) 

where a, b, and c are functions of temperature and we have specialized to SU{3). 
Taking 6 > 0, we see that if c were zero, the minimum of the potential would increase 
continuously from L = 0toL7^0asa went through zero. However, with the c term 
present \L\ will always iump discontinuously away from zero as a is lowered. In fact, 
no cases are known in which there is a second order phase transition involving an 
order parameter whose symmetries allow a cubic invariant in the free energy. The 
cubic term which is present because the symmetry is Z{3) forces the phase transition 
to be first order. 

The QCD phase transition is an inherently non-perturbative phenomenon. Per- 
haps the best long term strategy for studying it in full (non-universal) detail from 
first principles is to use lattice gauge theory. Following Wilson,lii the theory can be 
formulated on a lattice and the Euclidian path integral (|l|) can be evaluated using 
Monte Carlo techniques.0 This involves evaluating expectation values of operators by 
importance sampling — one chooses configurations weighted by a known probability 
measure — and is therefore limited to studying equilibrium situations in which no 
time evolution occurs. The prediction of Svetitsky and Yaffe described above that the 
pure SU (3) phase transition is first order has been confirmed in lattice gauge theory 



simulations. I 
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Including dynamical quarks in the partition function changes the situation en- 
tirely. (L) is nonzero even in the confined phase. This can be understood heuristically 
by noting that a static quark source can always be screened by virtual qq pairs in the 
ensemble. It can also be understood more formally^ by showing that after integrating 
out the fermions, the effective action for the gauge fields includes a term which breaks 
the Z{N) symmetry. No order parameter is known for the confinement /deconfinement 
transition in SU{N) gauge theory when dynamical fermions are present. This would 
seem to preclude the existence of a second order transition. A first order transition 
is logically possible, and presumably is obtained for heavy quarks. For two species of 
light quarks, on the other hand, lattice gauge theory simulations§@'§iii0iii'il 
show no signs of a first order deconfinement transition. L increases smoothly as a 
function of increasing temperature and no discontinuities are apparent. As in the 
ionization of a gas, there is no sharply defined transition, but rather a smooth change 
from one regime to another. 

Arguing only from fundamental symmetries, we have found that the phase tran- 
sition is first order in pure gluon QCD. We found an order parameter, but its sym- 
metry precluded a second order transition. In the presence of dynamical quarks, we 



were unable even to find an order parameter. Had we found a second order transi- 
tion, we could do much more. Powerful renormalization group techniques could be 
brought to bear on the problem, superceding the mean field analysis above. When 
there are massless quarks present, while there is no order parameter for the confine- 
ment / deconfinement transition, there is an order parameter for the chiral phase tran- 
sition. In Chapter 2, following arguments first advanced by Pisarski and Wilczekli 
which are similar in philosophy to those above, we will discover that for two species 
of massless quarks the chiral phase transition can be second order. This will have 
many implications and quantitative consequences. 

All the results of Chapter 2 describe quantities observable in thermal equilibrium 
where there is no time evolution. As we have discussed, this makes them ideal for 
comparison with lattice "experiments" . Working on the lattice, one has more freedom 
to adjust the parameters of the theory than in the real world. "Experimentalists" can 
choose what quark masses to use, while experimentalists cannot. We will see that 
this is unfortunate, because in order to observe the critical phenomena discussed in 
Chapter 2 on the lattice or in the real world, it will be desirable to use smaller up and 
down quark masses than nature does. This is possible on the lattice, although difficult 
in practice in current simulations. The future for "experimentalists" is bright. 

What about experimentalists? Experimentalists hope to study the physics of the 
QCD phase transition in relativistic heavy ion collisions. We will argue in Chapter 2 
that if the plasma in these events stays in equilibrium, it is unlikely that the chiral 
transition will have a dramatic signature. Of course, it is still possible to hope to 
detect the difference between T ^ T^. and T <C Tc. It is just that the physics at 
T Tc will be smooth. In studying heavy ion collisions, however, we really should be 
open to the possibility that the plasma cools through the transition without staying 
in thermal equilibrium. This will be our concern in Chapter 3. We will find that it 
is possible that long wavelength oscillations of the pion field are excited, leading to 
the production of clusters of pions within which the number ratio of neutral pions is 
far from 1/3. This is a potentially dramatic signature of the chiral phase transition 
in relativistic heavy ion collisions in which thermal equilibrium is not maintained. 

As previously mentioned, the second goal of this introduction is to sketch some 
standard ideas about what may happen in ultra-relativistic heavy ion collisions. We 
will not review this enormous field comprehensively. References to reviews wherein 
the interested reader can find both discussions of greater depth and further references 
will be given. The reader should also consult other articles in this volume. 

Let us recall a few qualitative facts which are known empirically from studying 
high energy proton-proton collisions with multiple production of secondary particles,B 
and which are presumed to apply to high enough energy collisions between nuclei 
also.iiiii The momentum distribution of the secondary particles which could in gen- 
eral be a function of the longitudinal and transverse momenta of the particles {pl,Pt) 
and the centre of mass energy of the incident protons {^/s) is in fact observed to be 
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independent of ^/s for ^/s larger than about 10 GeV and the invariant cross-section 
factorizes according to 

E, p, and m are the energy, momentum and mass of a particular secondary particle. 
The Feynman x variable is given by 

x = ^^^. (14) 

PL max V S 

The transverse momentum distribution F2 falls rapidly with increasing px and the 
mean transverse momentum is pt ~ 0.35 GeV.0 (Events in which hard parton-parton 
scatterings give rise to jets with large pt occur about once per million collisions.0 
Most high energy proton-proton collisions yield secondary particles with small trans- 
verse momenta.) To discuss the longitudinal momentum distribution it is useful to 
define a variable y, called rapidity, where 

1, ( E+pl \ , (E±Pl\ ,,,, 

One reason why y is convenient is that boosting in the beam direction corresponds 
simply to adding a constant to y. Another reason is that for p ^ m, 

?/ ~ ?7 = — logtan6'/2 (16) 

where the pseudorapidity rj is easily measured since it depends only on the polar angle 
6. A third useful property of y is that 

d!/=%. (17) 

We will work in the centre of mass frame, where pi = x = corresponds to y = 
and PL = ±Pl max corresponds to y = ±ymi,^ where 



1 

2 ~^\pt + 



ymax = -l0g( , . (1^ 



Feynman postulatecS that at large s there is a central rapidity region {\y\ <^ ?/max, 
or, equivalently, \pl\ <^ Pl max) where Fi{x) is a constant. In this region, 

EdV = Fi F2{pt) d\pT)dpL (19) 

and, integrating over pt and using (0), 

da 

— = constant . (20) 
dy 



Thus, there is a central rapidity region in which the multiphcity per unit rapidity 
is approximately constant. This result has indeed been seen in proton- (anti)proton 
collisions.0 The qualitative features found in high energy hadron-hadron collisions 
— low pt and a central rapidity plateau — are expected^ to characterize ultra- 
relativistic nucleus-nucleus collisions also. 

In a relativistic heavy ion collision of sufficiently high energy, the decay products 
divide naturally into two regions. In the centre of mass reference frame, the fragments 
of the incident nuclei are found at high (positive and negative) rapidity, while the 
central rapidity region is characterized by approximately constant multiplicity per 
unit rapidity and by low baryon number density.0 At early times, this central rapidity 
region experiences boost invariant longitudinal expansionill in which the momentum 
and position of particles are related in such a way that the rapidity is given by 



where t is the time and z is the position coordinate in the beam direction. Mean 
physical quantities depend only on the proper time 



and transverse position, and are independent of y. 

Describing the fragmentation regions requires consideration of QCD at nonzero 
temperature and baryon number density, and is not discussed in this article. We limit 
ourselves to addressing the physics of the "hot vacuum" in the central rapidity region. 
How high an energy is sufficient for there to be a reasonably clean separation between 
distinct fragmentation and central rapidity regions? The parameter of interest is 5y, 
the extent in rapidity of one of the fragmentation regions. 5y is not well known 
— estimates range from 4 to 8.11 Even if 5y ~ 4, present fixed target heavy ion 
collision experiments at the AGS and SPS do not access sufficient rapidities for there 
to be a baryon poor central rapidity region. At the AGS, ?/max ~ 2 and at the SPS 
Z/max ~ 3. In the coming generation of colliding beam heavy ion accelerators, the 
situation will be much improved. At the Relativistic Heavy Ion Collider currently 
under construction at Brookhaven National Laboratory, it is expected that Au-Au 
collisions with y/s = 200 GeV/A with multiplicities of at least 1400 per unit rapidity 
will be achieved.ii When the LHC is built at CERN and is run as a Pb-Pb collider, 
collisions with ^/s = 6300 GeV/A and with multiplicities of at least 2600 per unit 
rapidity will be possible.!! These multiplicity estimates are qualitative — they are 
based on extrapolating from data from proton-nucleus collisions, and results from 
event generators suggest that they could be too low.0 At RHIC (?/max ~ 5) the 
central rapidity region will have essentially zero baryon number density if Sy < 5 
and at the LHC (ymax ~ 9) there will be such a region for 6y < 9.11 Even if there 
are baryons in the central region at RHIC, the number ratio of mesons to baryons is 




(21) 



r = 



(22) 



8 



estimated to be ~ 10^.El (At the LHC, this ratio is expected to be > IO^.Ej) Thus, 
analyzing the physics of QCD at finite temperature and zero baryon number density 
will be of relevance in describing the physics of the central rapidity regions of heavy 
ion collisions at RHIC and the LHC. 

What happens in the central rapidity region? The standard description involves 
three phases — a pre-equilibrium phase, a hydrodynamic evolution phase during 
which thermal equilibrium is assumed, and a hadronization phase. There have been a 
number of different approaches to treating the earliest phase of a heavy ion collision. 
For example, in the string model, one imagines that the incident nuclei pass through 
each other and colour flux tubes, or strings, are drawn out between the receding frag- 
mented nuclei. These strings then break apart by quark-anti-quark pair production, 
and realizations of the scenario are based on Monte Carlo simulation of this process. 
This scenario has been reviewed by Bialas and Czyz.S It runs into difficulties at 
energies high enough that the string density becomes large. At high enough ener- 
gies, it should be possible to describe the earliest phase of the evolution in terms of 
parton interactions by perturbative QCD. This approach, called the parton cascade 
model, has begun to yield results. One starts with distribution functions for the 
partons in the nucleons of the incident nuclei, and follows the quark-gluon interactions 
in the framework of renormalization group improved perturbative QCD by solving 
transport equations for the parton distributions in phase-space with Monte Carlo 
methods. The results of this work give support at the parton level to the idea that 
a thermalized state is formed. Quarks thermalize more slowly than gluons because 
the relevant QCD cross-sections are larger by a factor of two or three.i3 Also, for 
both quarks and gluons, a thermal distribution of momenta is achieved more quickly 
than chemical equilibrium (equilibrium number density) is achieved. (It is not clear 
whether there will be enough time in a RHIC collision for quark chemical equilib- 
rium to be established. It is of interest that at the latest times considered by 
GeigeiH ( ~ 3 fm/c), the baryon number is in fragmentation regions of high rapidity, 
and the central rapidity region contains partons with thermalized momenta which 
are expanding longitudinally with little transverse motion. While these calculations 
are far from the end of the story of the pre-equilibration physics, they do provide 
justification for supposing that in an energetic enough colhsion (perhaps at RHIC) 
there is a low baryon number central rapidity region in which particle momenta are 
in local thermal equilibrium at T > Tc. 

If equilibrium is achieved, and while it is maintained, the evolution of the plasma 
can be described hydrodynamically. The matter is treated as a relativistic fiuid 
described by smoothly varying four-velocity, pressure, temperature, and energy den- 
sity without reference to the microscopic description in terms of partons. This ap- 
proach has a long history, but perhaps the first quantitative treatment was that of 
Bjorken.Eil The review by Blaizot and OllitraultH is a useful introduction to the sub- 
ject. Bjorken0 used the observation that the central rapidity region is approximately 
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boost invariant to find solutions to the hydro dynamic equations describing the lon- 
gitudinal expansion and consequent cooling of the plasma. Blaizot and Ollitraultil 
review various numerical simulations of the hydrodynamic equations which include 
transverse expansion. 

In using the hydrodynamical approach one assumes that thermal equilibrium is 
valid throughout the transition and until the hadrons no longer interact. "This is an 
assumption which is hard to justify in any satisfying way."i^ It is particularly difficult 
to see how long wavelength degrees of freedom can stay in equilibrium as the plasma 
expands and cools through the phase transition. Given the uncertainties, in consid- 
ering the physics of the chiral phase transition in relativist ic heavy ion collisions we 
will discuss both the phenomena that would occur if the plasma stays in equilibrium 
and those that could characterize a far from equilibrium transition. 

In the final phase of a heavy ion collision, regardless of what has happened earlier, 
hadrons, photons, and leptons impinge upon an experimentalist's detector. How is 
the experimentalist supposed to unravel what has gone on before? The zeroth order 
question is how can she determine whether for a brief instant a small region of quark- 
gluon plasma was created? Many signatures have been proposed for the production 
of a region of quark-gluon plasma in a heavy ion collision and its consequent passage 
through the QCD phase transition. All have problems to overcome, arising either from 
background, or from insufficient theoretical understanding, or from the possibility 
that the signature can be mimicked by hadronic processes. This suggests that the 
QCD phase transition will be detected by a process where various pieces of evidence 
become more and more convincing rather than by the observation of a "smoking gun" 
signature. If the speculations of Chapter 3 are correct, however, the situation could be 
improved. We will see that if the chiral phase transition occurs by a process that is far 
enough from thermal equilibrium that it can be idealized as a quench, then coherent 
long wavelength oscillations of the pion field may emerge. These would manifest 
themselves as clusters of pions with similar rapidity in which the number ratio of 
neutral pions is dramatically different from 1/3. If seen, this phenomenon would be 
a distinctive, qualitative signature that a far from equilibrium chiral transition had 
occurred. 



2. Critical Phenomena at a Second Order QCD Phase Transition 

The QCD phase transition has been studied from many different perspectives using 
many different models and approximations. Here, however, we will attempt to learn 
as much as can be learned analytically about the phase transition, relying as much as 
possible only on fundamental symmetries and universality arguments and as little as 
possible on specific assumptions and models. This philosophy yields little fruit when 
applied to the confinement / deconfinement transition. We saw in the previous chapter 
that in QCD with light quarks, this transition is neither first order nor second order 
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— it is a smooth crossover like the ionization of a gas. The chiral phase transition is 
another story. Wilczek0 has emphasized that in the chiral limit where there are two 
species of quarks with zero current algebra mass, the order parameter for the chiral 
phase transition has the same symmetry as the magnetization of a four component 
Heisenberg magnet, which has a second order phase transition. At an equilibrium 
second order phase transition, thermodynamic quantities have non-analytic behaviour 
determined by the physics of the modes which develop long correlation lengths at the 
critical point. Thus, by studying these universal characteristics in a much simpler 
model — the = 4 Heisenberg magnet — it is possible to learn much about the 
chiral phase transition in QCD. 

In this chapter, which to a large extent follows work done with Wilczek,! we 
explore the consequences of this approach to the QCD phase transition in thermal 
equilibrium. In the following section, we use renormalization group arguments to 
establish a dictionary between QCD and the magnetic system. Using this dictionary, 
we apply results for the static critical exponents obtained in the magnet model to 
QCD. In section 2.2, we discuss the behaviour of the pion and sigma masses at 
the transition. In section 2.3, we review recent work in lattice gauge theory which 
has begun to test these ideas. In section 2.4, we consider how the strange quark 
affects the phase transition, and again compare with lattice results. In the concluding 
section of this chapter, we discuss the implications of all this for cosmology and lattice 
gauge theory, and for heavy ion collisions under the assumption that the chiral order 
parameter remains close to thermal equilibrium through the phase transition. 

2.1. QCD and the 0(4) Magnet 

The physics of the QCD phase transition is qualitatively different in the cases 
of zero, one, two, or three or more flavours of quarks.0 In this section we consider 
QCD with two species of quarks. Pisarski and Wilczek^ used an analysis similar to 
the one which follows to show that for three or more flavours of massless quarks, the 
chiral phase transition is first order. In section 2.4 we will discuss the effects of a third 
quark species, massless or massive, and for now, therefore, we postpone our discussion 
of the effects of the strange quark. If there are two flavours of massless quarks, 
the QCD lagrangian is symmetric under global chiral transformations in the group 
SU (2) L X SU {2)r X U{1)l+r of independent special unitary transformations of the left 
and right handed quark fields, and a vector U{1) transformation which corresponds 
to baryon number symmetry. (The axial f/(l) which would make the symmetry 
group into U (2) x U{2) is a symmetry of the classical theory, but not of the quantum 
theory.0) This chiral symmetry breaks spontaneously down to SU{2)l+r x U{1)l+r 
at low temperatures, and is restored at sufficiently high temperatures. The order 
parameter for this phase transition is the expectation value of the quark bilinear 

M] ^ {cf^qR,) (23) 
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which breaks the symmetry when it acquires a nonzero value below some critical 
temperature Tg. 

At a second order phase transition, the system is at an infrared fixed point of 
the renormalization group. Renormalization group transformations correspond to 
coarsening one's view of the world. Thus, at an infrared fixed point, physics is scale 
invariant. The order parameter fluctuates on all length scales and, in particular, on 
arbitrarily long length scales. The correlation length associated with these scale in- 
variant fluctuations is necessarily infinite, and the correlation function which describes 
them is a power law. Our goal is not to describe all of the physics occurring at the 
transition. Rather, we limit our attention to the nonanalyticities of thermodynamic 
quantities. The partition function Z = Tre"^/^ is the sum of analytic terms, and so 
in a system with finitely many degrees of freedom must be an analytic function of the 
temperature. Nonanalyticity can only arise in the infinite volume limit, and so by 
limiting the scope of our treatment in this way we concentrate on the physics of the 
long wavelength degrees of freedom. At an infrared fixed point, renormalization [i.e. 
coarsening one's view) does not change the physics of the long wavelength fiuctua- 
tions of the order parameter, which is therefore independent of all short wavelength 
physics. The nonanalytic critical phenomena are universal: they only depend on the 
modes with long wavelength fiuctuations (and these are determined by the symme- 
tries of the order parameter) and on the scaling behaviour near the fixed point of low 
dimension operators constructed from these modes. 

In order to describe a second-order transition quantitatively, we must find a 
tractable model in the same universality class. For the chiral order parameter (PB| ) 
the relevant symmetries are independent unitary transformations of the left- and 
right-handed quark fields, under which 

M U^MV . (24) 

These transformations generate an U{2)l x U{2)ji symmetry, which is not quite what 
is needed, since it includes the axial baryon number symmetry which is not present 
in QCD. This problem is solved^ by restricting Ai to multiples of unitary matrices 
with positive determinant, instead of general complex matrices. Matrices A4 in this 
restricted class remain in this restricted class under the transformation (^^ only if U 
and V have equal phases. Hence, the axial U{1) has indeed been removed. The 2x2 
matrices Ai can conveniently be parametrized using four real parameters {a, n) and 
the Pauli matrices as 

M = (7 + in ■ T . (25) 
In fact the order parameter can be written as a four-component vector 

r = (a,7f) (26) 

and the transformations (0) are simply 0(4) rotations in internal space. Hence, 
the order parameter appropriate for the chiral phase transition in QCD with two 
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flavours of massless quarks has the symmetries of the standard 0(4) invariant = 
4 Heisenberg magnet. For smaller number of components, this sort of model is a 
much-studied model for the critical behaviour of magnets, with the order parameter 
representing the magnetization of a ferromagnet or the staggered magnetization of 
an antiferro magnet. 

In order to discuss the universal critical phenomena, it is sufficient to retain those 
degrees of freedom which develop large correlation lengths at the critical point. These 
are just long wavelength fluctuations of the order parameter, which is small in magni- 
tude near the critical point and therefore fluctuates at little cost in energy. Thus, the 
most plausible starting point for analyzing the critical behaviour of a second-order 
chiral phase transition in QCD is the Landau- Ginzburg free energy 

(27) 

Here /i^ is the temperature dependent renormalized mass squared. Tc is the tempera- 
ture at which /i^ vanishes, /i^ is positive above the critical temperature and negative 
below it. The quartic coupling A depends smoothly on the temperature. In this sec- 
tion we consider the case in which A is positive; we will discuss what happens if A is 
not positive at Tc in section 2.4. We neglect higher dimension operators, because they 
are irrelevant at an infrared fixed point — their coefficients go to zero upon renor- 
malization toward the infrared." The symmetry breaking pattern we want is (Ai) oc 1 
(equivalently, (a) ^ 0; (tt) = 0) below the transition which is indeed what we find at 
the minimum of the potential for positive A. This model has been studied in depth 
for arbitrary and spatial dimension d, and the existence of an infrared stable fixed 
point of the renormalization group has been established.ii Hence, it is a model for a 
second order QCD chiral phase transition for two massless quarks. 

Before proceeding, let us briefiy outline the the physics of the theory given by 




(|27|) in thermal equilibrium above, below, and at the critical temperature at which 
the second order phase transition occurs. At temperatures above the direction 
in internal space in which is oriented is correlated only over spatial distances of 
order a finite correlation length; over longer distances all orientations of are equally 
probable. Thus, is disordered and (0) = 0. At temperatures below Tc, is ordered: 
(0) is nonzero and selects a particular direction (defined to be the a direction) in 
internal space. Whereas the theory has an explicit 0(4) symmetry, the equilibrium 
configuration is not 0(4) symmetric — the symmetry is spontaneously broken. Fluc- 
tuations of about (0) in the a and vr directions are described by differing correlation 
functions. The explicit 0(4) symmetry implies that (0) could choose any direction 
in internal space. Below Tc, therefore, fluctuations of in the vf directions (which 
correspond to fluctuations of the orientation of (0)) are massless Goldstone modes. 

° (f)^ is not irrelevant in 3-dimcnsions, but as long as A is positive at Tc, (j)^ can be neglected relative 
to 0^ since |(/)| is small near the transition. The (ff' term will be important in section 2.4. 



13 



At T = Tc, (0) = 0, the correlation length is infinite, and the fiuctuations of (f) are 
scale invariant. 

We have established that there is an infrared fixed point of the renormalization 
group in the same universality class as QCD with two massless quarks. The theories 
in this universality class live in an infinite dimensional parameter space spanned by /i, 
A, and the coefficients of an infinite number of irrelevant higher dimension operators. 
The second order fixed point has some basin of attraction in this space of theories. 
For the rest of this chapter, we will explore the hypothesis that QCD is indeed in 
the basin of attraction of the infrared fixed point. We will calculate much about the 
QCD phase transition based on this assumption. We will also be show that it is 
consistent with many phenomena found to date in lattice gauge theory simulations. 
We can never be sure, however, that upon renormalization toward the infrared QCD 
is in fact driven to this fixed point. It is always a logical possibility that, say, QCD 
may live in a region of theory space from which infrared renormalization drives the 
theory into a region of negative A, making the phase transition first order. Hence, 
the strongest statement that it is possible to make is that it is plausible that in QCD 
with two species of massless quarks the chiral phase transition is second order. Our 
strategy, then, is to establish quantitative consequences of this hypothesis and test 
them against lattice simulations of QCD. 

When the free energy is written in terms of a and vf it looks much like the 
original model of Gell-Mann and Levy@ with two changes: there are no nucleon 
fields and only three (spatial) dimensions. These two changes refiect an important 
distinction.^ We are only proposing (P7| ) as appropriate near the second order phase 
transition point. This is because it is only there that we can appeal to universality — 
the long- wavelength behaviour of the a and vf fields is determined by the infrared fixed 
point of the renormalization group, and microscopic considerations are irrelevant to it. 
In Euclidean field theory at finite temperature, the integral over uo of zero temperature 
field theory is replaced by a sum over Matsubara frequencies uOn given by 2n7iT for 
bosons and {2n + l)7rT for fermions with n an integer. Hence, one is left with a 
Euclidean theory in three spatial dimensions with massless fields from the n = 
terms in the boson sums and massive fields from the rest of the boson sums and the 
fermion sums. Hence, to discuss the massless modes of interest at the critical point, 
(|27|) is sufficient and we do not need to introduce nucleon fields or constituent quark 
fields. In arriving at an effective three dimensional theory of the long wavelength 
fiuctuations of the order parameter near T^, all other bosonic degrees of freedom and 
all fermionic degrees of freedom have been integrated out. 

We have motivated a very definite hypothesis for the nature of the phase transition 
for QCD with two species of massless quarks, namely that it is in the universality class 
of the = 4 Heisenberg magnet. This means that under renormalization toward the 
infrared, QCD with two massless quarks is driven to the infrared fixed point of (|27| ) 
and that therefore the nonanalytic behaviour of thermodynamic quantities near Tc is 
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the same in the two theories. This hypothesis has numerous consequences which are 
the subject of the rest of this chapter. In the remainder of this section, we review the 
resulting predictions for the critical exponents.^ 

First, we define the reduced temperature t = (T — Tc)/Tc. The exponents a, 
f3, 7, 7], and u describe the singular behaviour of the theory with strictly zero quark 
masses as t — > 0. For the specific heat one finds 

C(T) ~ |t|-" + less singular. (28) 

The behaviour of the order parameter defines (3. 

(101) ~ \tf for t < . (29) 

rj and u describe the behaviour of the correlation length ^ where 

Ga^ix) = (0(x)„0(O)/3) - {(I)a){(pl3) ^ '^"/^ |^ | (d-2+q) ^^Pi'l^l/O ^^1 |a;| ^ OO . 

(30) 

A is independent of \x\, but may depend on t. The correlation length diverges as 
\t\ 0, and the correlation length exponent u is defined by 

^-l^r''- (31) 

Above Tc, where the correlation lengths are equal in the sigma and pion channels, the 
susceptibility exponent 7 is defined by 

£x G„/3(x) ~ r^. (32) 

We will discuss the behaviour of the susceptibility below the transition in the following 
section. At the critical point, the correlation length is infinite and the correlation 
function is a power law. The exponent r) is defined through the behaviour of the 
Fourier transform of the correlation function at t = 0: 

Gap{k ^ 0) ~ A;~2+'' . (33) 

The last exponent, 6, is related to the behaviour of the system in a small magnetic 
field H which explicitly breaks the 0(4) symmetry. Let us first show that in a QCD 
context, H is proportional to a common quark mass rriu = rrid = rriq. This common 
mass term may be represented by a 2 x 2 matrix T> given by rriq times the identity 
matrix. We are now allowed to construct the free energy from invariants involving 
both D and A4. The lowest dimension term linear in V is just trA^^P = niqa, 
which in magnet language is simply the coupling of the magnetization to an external 
magnetic field H oc nig. In the presence of an external field, the order parameter is 
not zero at Tc. In fact, 

{mt = 0,H^O)r~.H'/' . (34) 
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Thus when the 0(4) symmetry is exphcitly broken by a small external magnetic field 
or equivalently by small quark masses, the phase transition is a smooth crossover. It 
is, however, a very special crossover about which much can be said, because as H ^ 
one approaches the second order fixed point. 

The six critical exponents defined above are related by four scaling relations.^ 
These are 

a = 2 — du 

(3 = ^-{d-2 + ^) 

7 = (2-r/)z/ 

S = ^i^. (35) 

We therefore need values for r] and u for the four component magnet in c? = 3. 
These were obtained in the remarkable work of Baker, Meiron and Nickel,0 who 



carried the perturbative expansion in A for the theory ([27| ) to seven-loop order, and 
used information about the behaviour of asymptotically large orders, and conformal 
mapping and Pade approximant techniques to obtain 

T] = .03 ±.01 

u = .73 ± .02 . (36) 
Using (pSj), the remaining exponents are 



a = -0.19 ±.06 

p = 0.38 ±.01 

7 = 1.44 ± .04 

6 = 4.82 ± .05 . (37) 

Since a is negative there is a cusp in the specific heat at Tc, rather than a divergence. 
There are other ways in which these exponents could be calculated. In d = 4, mean 
field theory suffices and one has r] = and u = 1/2. We are interested in d = 3, 
however. A standard approach^ is to work in 4— e dimensions, evaluate the exponents 
order by order in e, and then set e = 1. The e-expansion has not been pushed to high 
enough order to compete in accuracy with the work of Baker, Meiron, and Nickel. 
Another approach is to do a finite temperature lattice simulation of the theory (|27|) 
and measure the critical exponents. This has recently been done to high accuracy by 
Kanaya and Kaya,0 who find u = 0.748 ± 0.009 and 7 = 1.48 ± .02, in agreement 



with the results of Ref. 45 



To reiterate, these exponents and the other critical phenomena we discuss in 
subsequent sections are universal. Whether the QCD phase transition is studied 
experimentally, or on the lattice, or by instanton liquid techniques, or by using the 
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Nambu-Jona-Lasinio model, or by using the (four dimensional) linear or nonlinear 
sigma models, or by any other means, if the hypothesis that the chiral transition is a 
second order transition in the 0(4) universality class is correct the exponents should 
turn out to have the values we have described. Of course, using any technique except 
appealing to the simplest model in the same universality class as we have done would 
make it much harder to study the theory near the critical point, and hence much 
harder to obtain the critical exponents. 

2.2. The Equation of State and the Pion and Sigma Masses 

The expressions which define (3, 7 and 6 are actually special cases of a more 
general relationship between the magnetization and the magnetic field called the 
critical equation of state. The equation of state has been calculated to order by 
Brezin, Wallace and Wilson.@ In this section, we review the use of the equation of 
state to determine the behaviour of the pion and sigma masses near the critical point. 

First, we must define what we mean by the "mass" of the pion and sigma. We 
could choose either to define the mass as an inverse correlation length or as an inverse 
susceptibility. We choose the latter, which is conventional in the condensed matter 
literature. Specifically, we define 



where (j)o = a and (pi = tTj, i = 1, 2, 3. Defining "masses" as inverse correlation lengths 
would give different scaling behaviours for masses as functions of t and H. (For rj = 0, 
however, both choices would yield the same scaling behaviour, and in the theory of 
interest rj is small.) We shall see that with the conventional choices (|38|) and (|39D, 
the masses can be extracted conveniently from the equation of state. It is worth 
noting that the masses we have defined are related only to the behaviour of spatial 
correlation functions in the static (equilibrium) theory. They carry no dynamical 
information, and should not be confused with, say, pole masses in a 4-dimensional 
theory. Also, we will only be able to make universal statements about how the masses 
scale at the transition. Normalizing the magnitudes of the masses will require using 
some specific model, and hence will not be universal. 

The equation of state gives the magnetization as a function of t and H. For 
the rest of this chapter, we will write the order parameter as M, for magnetization, 
keeping in mind that M = (a) = Also, recall that H, the external field, is 

proportional to the quark mass rrig. In order to define the equation of state, we first 
define a shifted field a = a — (a) = a — M. Then the equation of state is simply the 
relation 





and 




(39) 



(a) = . 



(40) 
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Fig. 1. From Ref. The order parameter M as a function of reduced temperature t = (T — Tc)/Tc 
for magnetic fields H = 0, 0.002, 0.005, and 0.02. Like t, both M and H are dimensionless. They are 
obtained from their dimensionful counterparts by dividing by non-universal dimensionful constants 
defined in such a way that for t < Q and H = the order parameter is given by M = (— i)'^, and for 
t = it satisfies M = H^/K 



This relation has been expanded to order by Brezin, Wallace and Wilson.c 
result can be expressed conveniently in terms of the variables 

y = H/M^ and x = t/M^'^ 



The 



(41) 



as 



y = fix) . (42) 

The function f{x) is universal and was calculated to order by Brezin, Wallace, and 
Wilson.^ Their result is also quoted in Ref. [3]. The units in which H and M are 
measured are chosen so that /(O) = 1 corresponds to t = and /(—I) = corresponds 
to the t < 0, H = coexistence curve. Knowing f{x), we can calculate the value of 
the order parameter M for a given H and t using (^2l). The behaviour of the order 
parameter is illustrated in Fig. 1. This figure and the other ones in this section 
should be viewed as illustrations of qualitative behaviour rather than quantitative 
predictions because they are based on setting e = 1 in the O(e^) expression for /(x). 
The values for the critical exponents themselves which we quoted in the previous 
section are quantitative predictions, complete with error estimates, because they are 
based on the much more elaborate analysis of Baker, Meiron and Nickel.Ei 
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From the equation of state, we can deduce the behaviour of and at nonzero 
(but small) t and H . The masses are given by 



^ dH 

< = m <«> 

and 

ml = —. (44) 

The first relation follows directly from the definition (|38|) , and the second follows from 
(|39|) and from assuming that M || H ^ so that a small change ± H gives a small 
change 5M ± M with 5M/5H = M/H. Using the equation of state, we can rewrite 
(H) and (H) as 

ml = M^-^f{x) (45) 



and 



ml = M'-' 1^5 fix) - ^f'ix)j . (46) 

Hence, 5 can be determined by measuring the ratio ml /ml at t = 0. In general, from 
f{x) we can find the pion and sigma masses for any t and H. 

There are two interesting limits which we will consider explicitly. First, for t > 
and if — > which corresponds to x — > oo, we should find the full 0(4) symmetry, 
and hence should find that the pion and sigma masses are identical. For x oo, the 
function /(x) behaves as f{x) = cx'^. The constant c is given to 0(e) in Ref. [3]. 
Applying (^5]) and (|46|) , we find that 

ml = ml = cf for a: — i> oo , (47) 

consistent with the symmetry. 

We can also consider the limiting case of approaching the coexistence curve. This 
means taking t < and if — > 0, which implies x —>■ —1. In this limit, M tends 
to a nonzero constant, and so from (^) , we obtain ml oc if, a familiar result for 
Goldstone bosons. The behaviour of the pion mass is illustrated in Fig. 2. 

The result (^) may look peculiar to a particle physicist who is more familiar with 
the zero temperature result 

"J = ^. (48) 



Before considering the sigma mass, we therefore pause here to explain how ( |18| ) and 
( |4^ ) are related. We have seen that ~ H and that the order parameter (qq) ~ 
(cr) ~ M. At zero temperature, is defined in terms of the axial current by the 
relation 

(0M^(0)|7r^(g))=zAg^5"^. (49) 
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Fig. 2. From Ref. [|. ml as a function of t for iJ = 0, 0.002, 0.005, and 0.02. Since M and H are 
in dimensionless scaled units, so is mi. For t = 0, ml ^ For = and t > (and for 

large enough t for any H) ml ~ f . For t <0 and H ^ 0, ml ^ H. 



(50) 



(51) 



In the zero temperature linear sigma model, the axial current is given by 
which means that /^r defined in (^8|) is simply 

A = (0 k 1 0) . 

This result suggests that we make the identification /^r ~ M, which does indeed make 
(|4^) and (48) equivalent. However, it is important to remember that using the linear 
sigma model at zero temperature cannot be justified by a universality argument in 
the way that using it near T = T^. can. Hence the argument of this paragraph is not 
a derivation of ( ^ID from the zero temperature result (^) . (|^) is valid near T = Tc 

0. Also, ttItt in (^8]) is a mass in a 3 -|- 1 dimensional Lorentz 
in (HI) is an inverse susceptibility in a 3 dimensional 



while (48) is valid at T 
invariant theory, while ml 
theory. We have simply shown that a reader familiar with one expression should not 
be surprised by the other. 

The behaviour of the sigma mass at the coexistence curve is trickier to obtain 
than that of the pion mass. First, we note that in mean field theory (e = 0) the 
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Fig. 3. From Ref. §. as a function of t ioi H = 0, 0.002, 0.005, and 0.02. Since M and H 
are in dimensionless scaled units, so is m^. For t — 0, = Srn^. The scale on the vertical axis is 
different than in Fig. 2. For H = and t > (and for large enough t for any H) = m^, and 
both increase like f . For t <Q and H ^ 0, the sigma mass decreases to zero as shown in Fig. 4. 

equation of state is simply y = f{x) = 1 + x, and m„ is easily evaluated using (^Gf ) . 
For if — > at fixed t < the result is 



Hence, in mean field theory decreases with H to a nonzero value at H = 0. 
However, for d < A when fluctuations are important, the result is quite different. In 
words, fluctuations of the massless pions produce new infrared singularities in the 
longitudinal susceptibility, or, equivalently, make the sigma massless. Now, let us see 
how this result can be obtained from the equation of state.Sii In the limit H —>■ 0, 
f{x) ~ H while f'{x), we will see, tends to zero more slowly. Hence, the second term 
in (^6]) is dominant and gives 



The difficulty is that for x — —1, f'{x) contains divergent terms like elog(x -|- 1), 
log^(x -|- 1) and log(x -|- 1). These terms do not exponentiate to f'{x) ~ (x -|- 1)^. 




(52) 




(53) 
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Fig. 4. From Ref. 1^. as a function of H ior t = —0.2. For large H it behaves as if it will 
be nonzero a,t H = 0, but in fact for _ff — > it decreases like ~ where to lowest order in e, 
1/2. 



Wallace and Ziap2l in fact find the result 

(j^f--^+-^y''- (54) 

Both the terms on the right side of ([5^) must be kept because they differ in their 
exponents only by order e. Also for this reason, the constants Ci and C2 calculated by 
Wallace and Ziat3 and given in Ref. |^ are only known to order e even though f{x) 
is known to order e^. 

Qualitatively, as H is lowered at fixed t < 0, at first the Ci term dominates and 
appears to be decreasing toward a nonzero value at = as in the mean field result. 
Then, the C2 term takes over and one finds that in fact the sigma mass goes to zero 
like oc H'^^'^. Recently, Anishetty et al^ have done an extensive analysis of the 
divergence of the sigma susceptibility produced by the massless pions for t < and 
H ^ 0. In agreement with the e-expansion result, they find that in d = 3, cx if^/^. 
The behaviour of the sigma meson mass is illustrated in Fig. 3 and Fig. 4. In future 
lattice simulations, as rriq is lowered toward zero, this behaviour should be observed. 
This result is an example of the power of the renormalization group techniques in 
obtaining universal results. If we had chosen a specific microscopic model, say that of 
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Gocksch,!! or the Nambu-Jona-Lasinio model of Hatsuda and Kunihiro,!! we would 
have been able to calculate non-universal quantities far from Tc, but would basically 
have been limited to using mean field theory, as those authors do. Then, we would 
have reached the incorrect conclusion that 7^ in the chiral limit below Tc. Here, 
by restricting ourselves to calculating universal quantities, we are limited to the region 
near the critical point, but our results are model independent and include the effect 
of fluctuations. 



2.3. Comparison with Lattice Simulations 

Finite temperature lattice QCD simulations are ideally suited to testing many of 
the predictions made in this chapter. For a pedagogical introduction to and review 
of finite temperature lattice QCD simulations, see DeTar's article in this volume. 
The static correlation functions of the three dimensional theory are natural objects 
to consider in finite temperature (Euclidean) simulations. Also, it is much easier to 
vary parameters like the temperature and the bare quark mass in a lattice simulation 
than in a real experiment. Hence, it should be possible to measure the static critical 
exponents of section 2.1, and the equation of state and the scaling behaviour of the 
pion and sigma masses of section 2.2 on the lattice. 

Present simulations Hffl'i'i'Sliiii provide evidence that in two flavour QCD 
the phase transition is second order, and that the order parameter is indeed A4. For 
example, in the results of Bernard et a/.@ shown in Fig. 5, there are no signs of 
any discontinuities in the expectation value of the order parameter as a function of 
temperature, (qq) decreases smoothly as a function of increasing temperature, as Fig. 
1 leads us to expect for simulations done with nonzero quark masses if the transition 
is indeed second order for niq = 0. All these authors have looked unsuccessfully for 
signals of two phases coexisting at one temperature, which would be characteristic 
of a first order transition. These simulations are consistent with the transition being 
second order, but do not rule out the possibility that in better simulations with smaller 
quark masses a small discontinuity could be found, indicating that the transition is 
in fact weakly first order. Using the zero temperature p mass to set the energy scale, 
Bernard et al^^ find that the temperature at which the decrease in the chiral order 
parameter occurs is about 140 — 160 MeV, significantly lower than in simulations 
without dynamical fermions. Over the same range of temperatures at which the 
chiral order parameter decreases, the Polyakov loop expectation value increases and 
the specific heat is large. Also, the calculations slow down — simulations must be 
run for a long time to get reliable results because large fluctuations occur. This is as 
would be expected near a second order critical point. Thus, in QCD with 2 dynamical 
quarks, it seems that the crossover associated with deconfinement, for which there 
is no order parameter, occurs at the same temperature as the chiral transition. The 
chiral transition is also a smooth crossover, but we expect it to approach a second 
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Fig. 5. From Ref. ||20| . The chiral order parameter in the lattice gauge theory simulations of Bernard 
et al. with two flavours of staggered quarks. The chiral order parameter is plotted as a function of T 
for niq — 0.025 and niq = 0.0125 in lattice units. The simulations were done on a 12^^ x 6 Euclidian 
lattice. On such a lattice, the period in the fourth direction is 1/T. Since the number of lattice 
sites in this direction (six) is not varied, to vary the temperature one varies the physical value of the 
lattice spacing. The fact that QCD is asymptotically free means that reducing the lattice coupling 
constant g reduces the physical value of the lattice spacing and therefore corresponds to increasing 
the temperature. Hence, in the figure temperature increases to the right. The curve labeled niq = 
is a linear extrapolation from the other two curves. We predict that if were lowered to zero, one 
would find the behaviour of Fig. 1 rather than a linear extrapolation. 



order transition as rriq is lowered further. All of this is in marked contrast to results 
of lattice gauge theory simulations with three or four species of light quarks, in which 
there is a first order chiral phase transition with latent heat, hysteresis, and abrupt 
changes in the values of observables.0'0'@'0 (Our results on the order of the chiral 
transition — second for two massless quarks and first for three or more — are also in 
agreement with results obtained by treating the QCD vacuum as an instanton liquid 
and making a mean-field approximation.^) 

Fig. 6. shows the behaviour of the screening mass (inverse correlation length) for 
the TT and cr in the same simulation as Fig. 5. If the exponent t] were zero, then 
and mo- (which, recall, are susceptibilities" ^Z^) would have the same scaling behaviour 
as the screening masses of Fig. 6. Since rj is in fact small, the screening masses can 
be viewed as crude approximations to and m^. As we expect, the sigma screening 
mass decreases and the pion screening mass increases as the temperature is increased, 
and above Tc the two masses each increase and appear to be becoming degenerate. 
This is evidence that the order parameter is indeed A4. 

So far, we have discussed various qualitative lattice results, that are in agreement 
with our expectations. Before turning to recent progress on more quantitative tests, 
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Fig. 6. From Ref. |^^. The inverse of the pion and sigma correlation lengths as a function of T 
for TOg = 0.0125 in lattice units in the same simulations as in Fig. 5. The screening masses are 
crude approximations to the inverse susceptibilities rn,r and Too- because rj is small. The qualitative 
behaviour should be compared to Figs. 2 and 3. 

it is worth enumerating the logical possibilities which could lead to a failure of the 
hypothesis that a second order chiral phase transition with 0(4) exponents will be 
seen in lattice simulations. Perhaps the most likely alternative is that the transition 
could in fact be first order. As we will see in the next section, this occurs if the strange 
quark is too light. More generally, we only know that an infrared fixed point with 
appropriate symmetry exists. We do not know that QCD is in fact in the basin of 
attraction of this fixed point. A second logical possibility is that the Ginzburg region 
may be too small to see the true critical behaviour. Fluctuations of the order param- 
eter are important only close to — for t = [T — Tc) /T^ < tc where to-, like T^, is not 
universal. Outside the Ginzburg region, [i.e. for t > to,) mean field theory is valid. 
Thus, even if the chiral phase transition does have 0(4) exponents, these exponents 
can only be measured by simulations done for temperatures satisfying t < to- Simu- 
lations with larger t would see the mean field exponents rj = and u = 1/2. to, like 
Tc, cannot be predicted by arguments based on universality, and must be measured 
on the lattice. Can to be very small? This is a logical possibility. For example]^ 
in ordinary, low temperature, BCS superconductors, to ~ {Tc/sfY where is the 
Fermi energy. In these materials, T^jef ~ 10~^ — 10~^ and Iq is minuscule. Recent 
numerical simulations by Kocc and K^gutB of the 2 + 1 dimensional Gross-Neveu 
model in which mean field exponents are seen can be interpreted as a sign that tc is 
small in this theory. In QCD, there is no analogue of the small parameter Tc/sf, and 
so there is no reason to expect to to be particularly small. Nevertheless, it remains 
a logical possibility. To exclude this scenario, it is therefore crucial to measure the 
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Fig. 7. From Ref. [E^. DeTar's plot of y vs. i using 0(4) exponents. 



critical exponents in QCD simulations with two flavours accurately enough to demon- 
strate that they deviate from their mean field values. The third logical possibility is 
that somewhere in theory space there is another infrared fixed point to which QCD is 
driven under renormalization. Were this the case, the chiral phase transition would 
be second order but would not have 0(4) exponents. Of the three logical possibilities, 
the third would be the most surprising. 

Recently, a number of authors have begun to attempt quantitative comparisons 
between lattice simulations and expectations based on the hypothesis that the chiral 
transition is second order with 0(4) exponents. DeTar has pointed out0 that it is 
possible to check that the order parameter M, the temperature t, and the symmetry 
breaking, term H ~ ruq are in fact related by a single equation of state (^). Following 
Binder ,Ea DeTar writes ( ^2] ) in a different, but equivalent, fashion. Define 

X = t/H^/^^ , y = M/H^'^ . (55) 

Then, the equation of state (^2]) is equivalent to 

y = gix) (56) 

where g{x) is a universal function related to f{x) of (^2|) by 

S: = g-\y)=x-'/^'f-\x) . (57) 
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Fig. 8. From Ref. |Q. Pseudocritical coupling as a function of quark mass for simulations done 
with N-r — 4 lattice sites in the Euclidean time direction. The dashed curve is a three parameter fit 
to (p9|). The solid curve is a two parameter fit with 1/(35 set to its 0(4) value. 

Fig. 7 is DeTar's plot of y versus x for a number of simulationS§S'i'§ with various 
quark masses, temperatures, and lattice sizes. If ( p6D is correct, all the points should 
lie on a single curve. In producing the figure, DeTar has used the 0(4) values for /3 
and (5, and has adjusted Tc to 150 MeV to obtain the best agreement, although he 
notes that 160 MeV works comparably well. DeTar's plot is encouraging, but not yet 
quantitative. For one thing, the preliminary data on the largest lattices with 12 lattice 
sites in the t direction are not in good agreement with the other data, suggesting that 
finite size effects may be important. In the future, with better simulations, it will be 
possible to apply this analysis more quantitatively — the simulation results will be 
fit to ( [55| ) vielding best fit values for Tc, (3 and 5, along with error estimates. 

KarschE^ has pioneered another way of measuring the critical exponents in a lattice 
simulation. From Fig. 3. we see that for H ^ 0, rricr has a minimum at t = tpc > 0, 
which Karsch calls the pseudocritical temperature. Because t, H, and M are related 
by (|4^) or equivalently (^), 

tpc ~ H'/^' . (58) 

Karsch therefore notes that in lattice simulations, the pseudocritical coupling Qpc at 
which the a susceptibility peaks should depend on rriq according to 

= Co + c,m'/^' , (59) 

where co is related to Tc for the second order transition with massless quarks, and 
neither co nor Ci are universal. Karsch and Laermann's plot0 of the pseudocritical 
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Fig. 9. From Rcf. @. The ratio of susceptibilities A = m^/m^ as a function of temperature for 
two different quark masses. 



coupling in various simulationsSSSi^ is shown in Fig. 8. Fitting the three param- 
eters Co, ci, and l/[35 yields l/[35 = 0.77 ± 0.14, which is slightly higher than the 
0(4) value = 0.55 ± .02 and which is in agreement with the mean field value 
l/PS = 2/3. However, fitting cq and ci with l/P6 fixed to its 0(4) value yields just as 
good a fit. Simple inspection of the figure leads to the conclusion that it is too early 
to extract the value of the exponent convincingly — simulations with lighter quarks 
are needed. 

Karsch and Laermann0 have evaluated both the tt and the a susceptibilities. 
From these, they extract a direct measurement of the exponent S. From (^) and 
(|46|) , we see that the quantity m^/m?^ which Karsch and Laermann call A is given by 

1 ml xfjx) 

For t > 0, A — > 1 for rriq — > and for t < 0, A ^ for ruq — >■ 0. Most important, 
however, is that for any quark mass, A = 1/5 at t = x = — i.e. at T = T^. Curves 
of A as a function of Q/g"^ for different quark masses should cross at Q/g"^, and at 
that coupling they should have A = 1/S. Karsch and Laermann's results are shown 
in Fig. 9. They find 

0.21 < 1/(5 < 0.26 . (61) 

This is in agreement with the 0(4) value 1/6 = 0.208 ± .002 and disagrees with the 
mean field value 1/5 = 1/3. To make this result even more compelling, it would be 
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good to have curves for lighter quarks to see that they all cross at the same Tc and 
A, and to see A decreasing toward zero as rriq — > for T < T^. 

One critical exponent which has not been mentioned in this section is a, which is 
associated with the specific heat. Because a < 0, there should be a cusp in the specific 
heat. However, the specific heat can also have a large analytic contribution. In fact, 
as the temperature is raised through the chiral transition, deconfinement also occurs 
and the number of degrees of freedom in thermal equilibrium increases substantially, 
leading to a large specific heat.iil This large contribution to the specific heat is an 
analytic function of T even for = 0, while the contribution from the chiral order 
parameter is not. Nevertheless, because the smooth part of the specific heat is so 
large, extracting the non-analytic part and thus a will be difficult.E3 

By this point it should be clear that quantitative tests of the hypothesis that 
the QCD phase transition is in the same universality class as the 0(4) magnet are 
just beginning to be possible. To date, lattice simulations of two flavour QCD are 
in qualitative agreement with this hypothesis. As lattice simulations improve, these 
tests will become more and more quantitative. It should soon be possible to exclude 
conclusively the possibility that the exponents take on mean field values. Indeed, the 
results of Karsch and Laermann for 1/6 are close to doing this already. This would 
demonstrate that the simulations are close enough in temperature to Tc that the 
fluctuations of the order parameter are important — i.e. the simulations are within 
the Ginzburg region. This would be reassuring, as it would imply that as simulations 
improve further the true critical behaviour will be explored, and the hypothesis that 
the transition is in the 0(4) universality class will be tested quantitatively. 

Now, some caveats. Present lattices are small. Ideally, one wants the correlation 
length to be large compared to the lattice spacing and small compared to the lattice 
size, and this requires much larger lattices. (Berera^ has suggested, however, that 
this problem can be turned into a virtue — finite size effects and finite lattice spacing 
effects affect the critical behaviour and therefore can be used to measure critical 
exponents like u.) However, the fundamental reason why present simulations do 
not offer more quantitative tests of the predictions in this chapter is that to date 
quark masses have been so large that correlation lengths do not get very long at 
Tc. For example, in the work of Bernard et al.^ the correlation length in the pion 
channel at Tc is only about 2.5 lattice lengths. To really probe the critical phenomena, 
lighter quarks are needed. Unfortunately, longer correlation lengths are necessarily 
accompanied by numerical critical slowing, and this makes simulations challenging. 

Recent work by the Columbia group" offers another indication that lighter quark 
masses are necessary in order to quantitatively probe the critical phenomena. They 
vary the "valence quark mass" (the quark mass in the qq operator inserted in order to 
measure the order parameter) and the "sea quark mass" (the mass of the dynamical 
quarks in the fermion determinant) independently, rrisea cannot be made too small 
because to do so would slow the simulation down; mvai can be made arbitrarily small. 
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The first result they find is that the order parameter goes to zero as a temperature 
dependent power of m^ai for fixed mgea = 0.01. (All masses are in lattice units.) They 
interpret the power law behaviour they see for mvai between about 10~^ and 10~^ 
as a sign of some sort of critical behaviour, but no interpretation in terms of the 
exponents we have discussed is possible, since these exponents are only relevant for 
the physical case where mvai = msea, = frtq- They also find indications that for fixed 
mvai, variations of msea which are accessible with their present simulations {i.e. down 
to m-sea = 0.01) seem to renormalize Q/ and have no other effect. (It is knownil that 
the effect of heavy quarks is to shift 6/g'^.) If these preliminary results are confirmed, 
they suggest that mgea must be reduced further than has been possible to date in 
order to see the critical phenomena characteristic of — >■ discussed in this article. 

Another hurdle to be overcome before lattice simulations can measure the critical 
properties of the QCD phase transition is that any lattice implementation of fermion 
fields only exhibits the full chiral symmetry in the continuum limit. Thermodynamic 
simulations with Wilson fermions are difficult because Wilson fermions break chiral 
symmetry completely. Larger lattices or an improved action are required to avoid 
lattice art if acts. mi If staggered fermions are used, one starts with four fiavours of 
fermions which in the continuum have an SU (4) x SU{A) chiral symmetry. On the 
lattice, a U{1) x f/(l) subgroup is all that remains. In order to study two fiavours of 
fermions, one takes the square root of the fermion determinant in the lattice action. 
It is not at all clear what this does to the lattice chiral symmetries. What is known is 
that all three pions become light only in the continuum limit. Finite lattice spacing 
effects leave one pion light, while making two of them heavy. A transition in which a 
U{1) xU{l) symmetry breaks to U{1) would be in the universality class of the N = 2 
magnet, as Boyd et alM have discussed. It is therefore possible that 0(2) exponents 
will be measured at first, and the 0(4) exponents will be seen only as finer lattices 
become possible and the continuum limit is approached. Therefore, it would be nice 
to be able to distinguish 0(4) exponents from 0(2) exponents, but this will require 
much greater accuracy than is presently possible. (For example, 6 = 4.808 ± .007 and 
l/l36 = 0.60 ± .01 at an 0(2) transition, and whereas distinguishing 0(4) exponents 
from mean field exponents seems feasible in the short term, distinguishing between 
0(4) and 0(2) will be more difficult.) Therefore, for staggered as for Wilson fermions, 
before we are able to make confident quantitative tests of our predictions for the 
critical phenomena we need finer lattices so that extrapolation to the continuum 
limit can be done. Eventually, the results for staggered and Wilson fermions must 
agree when so extrapolated. 

At present, it is clear, many caveats remain to be dealt with before quantitative 
measurements of the critical exponents in lattice simulations become compelling. The 
situation will improve steadily as simulations become possible on larger lattices with 
smaller quark masses. In the future, simulations which have longer correlation lengths 
and which are closer to the continuum limit will be able to measure critical exponents. 



30 



correlation functions, susceptibilities, and the critical equation of state, and test the 
hypothesis that the chiral phase transition in QCD with two flavours of quarks is 
in the same universality class as the 0(4) magnet. At present, this hypothesis has 
passed a number of qualitative tests. 

2.4- The Influence of the Strange Quark 

To this point in this chapter, we have described a world with two massless or 
light quarks, and hence we have implicitly been taking the strange quark mass to 
be infinite. Pisarski and Wilczek showedca that if the up, down, and strange quarks 
are all massless, then the chiral phase transition is first order. With three massless 
quarks, the order parameter A^* is a 3 x 3 matrix and there are two possible quartic 
couplings: Aitr(7W''7Vl)^ and X2{tTA4'' Ai)"^ . The renormalization group equations for 
this model have been studied to lowest order in e by Pisarski and Wilczekii and 
PatersonJ^ There is indeed a fixed point, but it is not stable. For three or more 
flavours of massless quarks, there are always directions in Ai — A2 space around the 
fixed point for which renormalization toward the infrared drives the theory away from 
the fixed point. Assuming the 0(e) result that there is no infrared stable fixed point 
is correct, the phase transition must be first order. This result has been verified by 
Gausterer and Sanielevici.i^ They did numerical simulations of the 3x3 matrix model 
and found that the transition is indeed first order. 

The chiral phase transition is second order in a world with two massless quarks 
and first order in a world with three. Hence, as the strange quark mass is reduced 
from infinity to zero, at some point the phase transition must change from second 
order to first order. This point is called a tricritical point. In a lattice simulation, the 
strange quark mass could be tuned to just the right value to reach the tricritical point. 
We now discuss the critical exponents that would be observed in such a simulation. 

Let us consider the effect of adding a massive but not infinitely massive strange 
quark to the two flavour theory. This will not introduce any new fields which become 
massless at Tc, and so the arguments leading to the free energy (^) are still valid. The 
only effect of the strange quark, then, is to renormalize the couplings. Renormalizing 
fi"^ simply shifts Tc, as does renormalizing A unless A becomes negative. In that case, 
one can no longer truncate the Landau- Ginzburg free energy at fourth order. After 
adding a sixth order term, the free energy becomes 



While for positive A, 0^ increases continuously from zero as /i^ goes through zero, 
for negative A, 0^ jumps discontinuously from zero to |A|/(2k) when fi^ goes through 
A^/(4k). Hence, the phase transition has become first order. Thus, at some rris = m* 
at which A = at T = Tc, the phase transition changes continuously from second 
order to first order. 




(62) 
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The nonanalytic behaviour of thermodynamic functions near tricritical points, 
just as near ordinary critical points, is universal. Hence, it is natural to propose^y 
that QCD with two massless flavours of quarks and with T near Tc and near 
its tricritical value m* is in the universality class of the 0^ Landau-Ginzburg model 
This model has been studied extensively.0 Because the 0^ interaction is strictly 
renormalizable in three dimensions, this model is much simpler to analyze than the 
0"^ model of the ordinary critical point. No e expansion is necessary, and the critical 
exponents all take their mean field values. There are calculable logarithmic corrections 
to the scaling behaviour of thermodynamic functions,0 but we will limit ourselves 
here to determining the mean field tricritical exponents. 

In mean field theory, the correlation function in momentum space is simply 
Gapik) = dapik"^ + /i^)""*^- Since jj? ~ t, this gives the exponents ?7 = 0, 7 = 1 
and V = 1/2. To calculate a and /?, we minimize F for H = X = V0 = 0, and find 
a = 1/2 and /5 = 1/4. To calculate 6, we minimize F for t = X = V0 = and find 

5 = 5. 

The result for the specific heat exponent a is particularly interesting, since it 
means that the specific heat diverges at the tricritical point, unlike at the ordinary 
critical point. This means that whereas for large enough that the transition is 
second order the specific heat C(T) has a cusp but is finite at T = Tc, as is lowered 
to the tricritical value C{Tc) should increase since at the tricritical point it diverges. 
This behaviour should be seen in future lattice simulations. 

At a tricritical point there is one more relevant operator than at a critical point, 
since two physical quantities {t and m^) must be tuned to reach a tricritical point. 
Hence, a new exponent (pt, the crossover exponent, is required. For A 7^ 0, tricritical 
behaviour will be seen only for |t| > t*, while for |t| < t*, either ordinary critical 
behaviour or first order behaviour (depending on the sign of A) results, t* depends 
on A according to 

t* ~ A^/"^* (63) 

The mean field value of (pt is obtained by minimizing the free energy F for H = V0 = 
0, and is 0t = 1/2. These mean field tricritical exponents, a = l/2, /5 = l/4, 7 = 1, 

6 = 5, 7] = 0, h> = 1/2, and 0t = 1/2 describe the real world if is close enough 
to the tricritical value, and will describe future lattice simulations with chosen 
appropriately. 

We now construct a phase diagram showing the order of the chiral transition as 
a function of rUg and the light quark mass Trig. In the {mq,ms) plane, there is a 
region around (00, 00) in which all quarks are sufficiently massive that a first order 
deconfinement transition is obtained. In a region around (0,0), all quarks are light 
enough that the first order chiral transition characteristic of the theory with three 
massless quarks is obtained. On the left side of Fig. 10, nig = and the chiral phase 
transition is first order for rris < m* and second order for > m*. For rUg 7^ 0, 
and ms large enough, the chiral transition is a smooth crossover. For > m* 
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Fig. 10. A phase diagram showing how the order of the chiral transition depends on and iris. 
For large m^, there is a second order transition with 0(4) exponents at ruq = 0. For ruq nonzero, 
there is a smooth crossover which, for smah rriq, looks like one of the curves of Fig. 1. For niq = 0, 
the transition becomes first order when is reduced below its tricritical value m*. In a region 
around rUs = rriq = 0, the transition is first order. The boundary of this region is a curve of second 
order transitions with Ising model exponents.O In the figure, this curve is drawn with the shape 
(to* — TOg) ^ TOq^'^, which is valid near the tricritical point. 



and nig nonzero but small enough, the crossover is described by 0(4) exponents, as 
discussed in previous sections. At the boundary between the first order region and 
the smooth crossover region, the transition has been analyzed in mean field theory 
by Gavin, Gocksch and Pisarski.llll They find that along this curve in the {mq,ms) 
plane, a single scalar field is massless and therefore conjecture that along this curve 
the transition is second order with Ising critical exponents. (The scalar which is 
massless is pure ss near nis = 0, becomes an SU{3) singlet as one follows the curve 
to the point at which = nig and becomes the a of previous sections as one follows 
the curve further toward the tricritical point at which niq = 0.) Using (|62D, we 
can deduce the shape of the curve of second order transitions near the tricritical 
point. For A < 0, consider turning on a small H. H breaks the 0(4) symmetry, and 
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makes the order parameter point in the a direction. Working to lowest order in H, 
we find that a jumps discontinuously between degenerate minima at a = 2Hk,/\^ 



and o = J\X\/{2k) — {1/2)Hk/ \^ at the temperature at which jj? goes through 



A^/ (4k) + iJy 2/t/l A|. As H increases, the size of the jump in a at Tc shrinks, until at 
some H ~ |A|^/^/k^/^ the two degenerate minima merge and the transition is second 
order. For even larger H, the transition is a smooth crossover. Since H ~ iriq and 
A ~ {nis — m*), we conclude that near the tricritical point, the curve of second order 
transitions separating the region in which the transition is first order from the region 
in which it is a smooth crossover has the shape 



Away from the tricritical point, a curve of second order transitions exists,^ but 
nothing is known about its shape. 

The real world is one point on Fig. 10. Lattice simulations, on the other hand, 
can be done with varying quark masses. Thus, on the lattice unlike in the real world 
it will be possible to test and explore all the physics of Fig. 10. It is clearly very 
important to know where the physical values of rris and rriq lie in Fig. 10. This is not 
a universal question, and therefore can only be answered by experiment or by lattice 
simulations. 

In simulations done on a 16'^ x 4 lattice with staggered fermions,il the Columbia 
group found a first order transition for = = 0.025 in lattice units, and found no 
sign of first order behaviour in simulations with only two fiavours of fermions and in 
a simulation with = 0.1 and rriq = 0.025. By measuring zero temperature hadron 
masses, they estimate that in the latter simulation mq is larger than its physical value 
and, most important, is smaller than its physical value. These results suggest that 
in this simulation, rris is greater than m*, and furthermore that the physical value of 
rris is even larger. Thus, the Columbia results suggest that the real world is at a point 
in Fig. 10 somewhere in the region labelled "crossover described by 0(4) exponents." 

Recently, Iwasaki et a/.0 have announced results from simulations with two light 
and one heavy Wilson fermions. In agreement with the Columbia results, they find 
a first order transition when all three quarks are light, and no sign of a first order 
transition in simulations with two light quarks. However, they report seeing a two 
state signal characteristic of a first order transition in simulations with two light 
quarks and a strange quark either of mass ~ 150 MeV or of mass ~ 400 MeV. 
This suggests that the real world lies in the region of Fig. 10 labelled "first order" , 
and is inconsistent with the Columbia results. There are a number of extensions 
to the work of Iwasaki et al. that would make it more compelling. First, in order 
to demonstrate consistency with the result (from their own simulations and those of 
others) that the transition is a smooth crossover for two fiavours, it would be good if a 
simulation was performed with a heavy enough strange quark that no signs of a first 
order transition are seen. Second, in order to more confidently make comparisons 





(64) 
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with the Columbia results, a measurement of the zero temperature hadron masses 
would be helpful. Finally, as discussed briefly in the previous section and at length 
in DeTar's article in this volume, studying the chiral transition with Wilson fermions 
is tricky because of potential lattice artifacts since the chiral symmetry is completely 
broken on the lattice. At the end of the day, however, when simulations are done on 
large enough lattices, it should not matter whether staggered or Wilson fermions are 



used. Although the results of Ref. |^ are at present more convincing than those of 



Ref. ||72|, the question of where on Fig. 10 the real world lies remains open. Time 
(i.e. further simulations) will tell. 



2.5. Implications 

We have seen in the previous two sections that the physics of the chiral phase 
transition in QCD is well-suited for study on the lattice. The hypothesis that the 
chiral transition in two flavour QCD is second order with 0(4) critical behaviour 
has many implications which are being and will be tested in lattice simulations. The 
exploration of the (m^, m^) plane has also begun. Let us now turn to real experiments, 
as opposed to those on the lattice. Here, we are not free to dial the bare quark 
masses. We will see that this is unfortunate, particularly in our discussion of heavy 
ion collisions. First, however, let us dispose of two other possible arenas for testing 
our results. 

The QCD phase transition certainly occurred in the early universe. Indeed, much 
work has been done on possible observable cosmological effects of this transition, if it is 
first order. There have been speculationsEl that a strongly first order transition could 
have affected big bang nucleosynthesis and therefore left a signature in the primordial 
nuclear abundances,^ or could have left an imprint in the form of gravitational waves, 
or could even have left relic chunks of new forms of matter. Unfortunately, we have 
seen that for physical values of the strange quark mass, it is likely but not certain 
that the transition is second order. In fact, because the up and down quarks are 
not massless, the second order transition is smoothed into a crossover. We can think 
of no observable consequences of a second order QCD phase transition in the early 
universe, and even fewer consequences of a smooth crossover. 

There is one possibly useful cosmological (non) consequence of a QCD transition 
which is not first order. In the early universe, the expansion rate during the transition 
is very slow compared to QCD time scales, and thermal equilibrium is maintained. 
Therefore, if the transition is not first order there is no entropy production at the 

''A strong first order phase transition creates, inhomogeneities in the baryon number density of the 
universe which later affect nucleosynthesis.O Present observations of nuclear abundances actually 
constrain big bang nucleosynthesis tightly enough that they disallow models irLjiSfhich the abundances 
are significantly affected by inhomogeneities in the baryon number density.O These cosmological 
observations, therefore, are evidence that the QCD phase transition is not strongly first order, in 
agreement with our results. 
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transition. Thus, whatever mechanism generates the baryon number density of the 
universe at some higher energy scale need only generate the present baryon to entropy 
ratio, and need not generate a larger baryon number to compensate for dilution at a 
first order QCD phase transition. 

Now, let us turn from cosmology to laboratory experiments. Since certain anti- 
ferromagnets including dysprosium have order parameters in the = 4 universality 
class,Ei experiments on thephase transition in these materials could help us under- 
stand the QCD transition.eO Alas, in these magnets there is a quartic operator like 
{nf + 7r|)(7r3 + o"^) which is allowed by the microscopic hamiltonian of the magnets,@ 
but not by that of QCD. This operator makes the symmetric Heisenberg fixed point 
infrared unstable, and either makes the phase transition first order or makes it second 
order but governed by an anisotropic fixed point. Hence, dysprosium is not a suitable 
arena for learning about the QCD phase transition. 

Finally, we turn to relativistic heavy ion collisions. As we discussed in the Intro- 
duction, in a sufficiently high energy collision the baryon number of the incident inci- 
dent nuclei ends up in that part of the plasma heading approximately down the beam 
pipes, and the central rapidity region consists of a hot plasma with approximately 
zero baryon number. In the remainder of this article, we will study the behaviour of 
the plasma in the central rapidity region as it expands and cools through T = Tc and 
eventually hadronizes and becomes pions which fly off and are detected. 

The defining characteristic of a second order phase transition is the divergence of 
correlation lengths. How could this feature be observed here? Large volumes of space 
with the order parameter correlated and pointing in a direction different than the 
true vacuum {i.e. sigma) direction will become regions in which the order parameter 
oscillates coherently about the sigma direction. After hadronization, correlated vol- 
umes will turn into regions of space where the ratio of the number of charged pions to 
neutral pions has some fixed value. In the standard scenarioiil described in the Intro- 
duction, the longitudinal expansion is boost invariant and different positions in the 
plasma along the beam direction become different rapidities as the plasma expands. 
Therefore, one would hope that a signal of a second order phase transition would be 
fluctuations in the ratio of charged to neutral pions as a function of rapidity. Such 
phenomena will be our subject in the next chapter. 

A prerequisite for the fluctuations discussed above to be observable is that the 
correlation length must get long compared, say, to Tc. To determine whether this 
does indeed happen, we must leave our universality safety net behind, since neither 
Tc nor the magnitudes of correlation lengths are universal. From Fig 2, it is im- 
mediately apparent that we have a problem. The longest correlation length is in 
the pion channel, and the pion mass is increasing with temperature. This suggests 
that the pion mass at Tc is larger than irtT^iT = 0) = 135 Me V. This result can be 
checked by comparing with (non-universal) results obtained in a variety of models. 
In chiral perturbation theory using the nonlinear sigma model,0 in the linear sigma 
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model,E§ and also in the Nambu-Jona-Lasinio model,!! the pion mass increases from 
its zero temperature value as the temperature is increased from zero. Working in a 
specific four-dimensional model forces one to make a mean fieldi or one-loop0 or low 
temperature0 approximation which is not valid near T = Tc. This precludes obtain- 
ing good results for universal quantities, but allows rough estimation of non-universal 
quantities like and the low temperature behaviour of the pion mass. Non-universal 
analyses suggest that 777-7^ increases with temperature as T increases from zero, and 
our universal result shows that it is increasing with temperature for T ~ T^. Hence, it 
seems clear that the longest correlation length at will be shorter than (135 MeV)~^. 
This is to be compared to Tc itself, which for the case of two massless quarks is around 
140-160 MeV.Byc3 Hence, even though the quark masses are indeed small (~ 10 MeV), 
the magnetic field H proportional to iriq is large enough to prevent any correlation 
lengths from reaching interesting values. 

There is an appropriate quantitative criterion to determine whether an equilibrium 
second order phase transition leads to dramatic effects. One compares the energy in 
a correlation volume just below with the zero temperature pion mass to determine 
whether or not the correlated volume can become a large number of pions. Using 



the lattice simulations of Ref . pO[ , we can make a crude attempt at this comparison. 
The sum of the energy and pressure in a correlation volume is about 1/4 the zero 
temperature p mass. Taken literally, this means that each correlation volume becomes 
only one or two pions in the detector. In these simulations, iriq is somewhat larger than 
its physical value, perhaps by a factor of two, and there are many other caveats as we 
have seen, so this estimate should not be taken literally. Nevertheless, it seems clear 
that the physical value of niq is large enough that in an equilibrium phase transition 
a correlation volume at does not evolve into a large number of correlated zero 
temperature pions. This is not encouraging. 

We pause here for an aside. The reader may be wondering why, when the seem- 
ingly small equal quark mass m„ = = rriq has such deleterious effects, we have 
completely neglected the difference between the up and down quark masses. Unequal 
quark masses allow terms of the typeS 

AF oc (5m) 2 (a^ -A + '^1 + A) ■ (65) 

If one is close enough to the critical point that this term matters, one will discover an 
anisotropic fixed point rather than the symmetric Heisenberg fixed point. However, 
while the effect of a common quark mass, namely the mass of the pion, is comparable 
to Tc, the effect of (|65|) is much smaller. For example, the QCD contribution to the 
difference in mass be^een the charged and neutral kaon= i= about 5 MevS Therefore, 
we need not worry about (^5]) in real experiments since it is much less important than 
the effect of the "magnetic field" proportional to the common up and down quark 
mass. Of course, ( |65D could be introduced and studied on the lattice. 

It seems likely that if the standard scenario for heavy ion collisions in which 
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the QCD plasma cools through Tc while staying close to thermal equilibrium is cor- 
rect, then no correlation lengths will get long enough for there to be any dramatic 
observable effects of the phase transition. The chiral "transition", like the confine- 
ment /deconfinement "transition," will be a smooth crossover. If we were able to dial 
down the quark masses and hence the pion mass, phenomena associated with a second 
order transition would become more prominent. Alas, in the real world, unlike on a 
lattice, we have no such freedom. Let us hope that we will reach less disappointing 
conclusions in the next chapter when we consider the physics of the QCD phase tran- 
sition in a heavy ion collisions in which the chiral order parameter does not stay in 
thermal equilibrium. 

The future of the study of the chiral phase transition in QCD on the lattice 
looks promising. As simulations improve, investigations of the plethora of static 
critical phenomena we discussed in the first sections of this chapter will become better 
and better. Critical exponents, the equation of state, the critical behaviour of the 
pion and sigma susceptibilities, and tricritical exponents, are all out there waiting 
to be measured. The hypothesis that QCD with two flavours of massless quarks 
has a second order chiral phase transition in the same universality class as the 0(4) 
Heisenberg magnet will be tested more and more quantitatively. 



3. Prom Disorder to Long Wavelength Oscillations: The QCD Phase Tran- 
sition Par Prom Thermal Equilibrium 

The gloomy penultimate paragraph of the previous chapter began with a condi- 
tional sentence. In this chapter we will consider the observable effects in heavy ion 
collisions if the plasma does not stay close to thermal equilibrium through the tran- 
sition. There are tantalizing hints in cosmic ray physics that point in this direction. 

Among the zoo of high energy cosmic ray events which have been detected in 
emulsion experiments are a particularly peculiar class of events called CentauroJ^^ 
discovered in emulsion experiments at the Mount Chacaltaya Observatory. Centauros 
are events with total energy of order 1000 TeV in which many (of order 100) charged 
hadrons each with energies of a several TeV and very few photons or electrons are seen 
in a cosmic ray induced shower. In a normal event, the charged hadrons are accom- 
panied by photons from the decay of neutral pions. Centauros have also been seen in 
emulsions exposed at the Pamir Observatory and the Mount Fuji Observatory. § The 
experimental situation is murky, however, because more recent searches at Mts. Kan- 
bala and Fujil^ failed to find any Centauros.ii Events called mini- Centauros which 
have fewer charged hadrons than Centauros but which also have very few photons or 
electrons have also been observed. In the sample of events in described by Lattes et 
a/.,0 there were 5 Centauros and 13 mini-Centauros, representing about 1% and 2% 
of the events seen with energy of the appropriate order of magnitude. Cosmic rav 
shower simulations do predict events which can be interpreted as mini- Cent aurosj^ 
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but it is difficult to interpret Centauros normally.E3 

The JACEE collaboration has flown balloon borne emulsion detectors to study 
cosmic rays.ii They study cosmic ray interactions in which the primary scattering 
occurs within one of the many layers of the detector. This is a big advantage, because 
it means the location of the vertex is known to within a few microns, the tracks of 
all the decay products are well measured, and therefore the pseudorapidities of the 
particles are known. The disadvantage, however, is that balloon borne detectors are 
necessarily exposed for much shorter periods of time than detectors left on mountain 
tops. Mountain top experiments therefore detect higher energy (rarer) cosmic rays. 
In fact, JACEE has seen very few events with energies in the 1000 TeV range and 



no Centauros.E3 However, recently!^ they have discovered events in which there are 
"ant i- Cent auro" regions in pseudorapidity-azimuth phase space containing a few tens 
of photons and almost no charged tracks. 

Centauros are peculiar because so many charged pions are observed without any of 
the gammas that would indicate the decay of neutral pions. This apparent violation of 
isospin invariance is puzzling, unless one thinks of it in the language of a second order 
phase transition in which these events can be interpreted as the creation of a volume 
of QCD plasma in which the field has fluctuated throughout most of the plasma in 
some direction in the vri — 1^2 plane. This implies correlation throughout most of a 
volume of plasma large enough that it becomes about 100 zero temperature pions. We 
convinced ourselves above that this could not happen if the plasma remains close to 
thermal equilibrium. Hence, these Centauro events provide a tantalizing hint that it 
might be wise to consider the effects of going from the disordered phase above to the 
ordered phase rapidly without maintaining thermal equilibrium. Regardless of what 
the final word on Centauro cosmic ray events turns out to be, it is worth considering 
qualitative phenomena which can arise if the chiral transition is far from thermal 
equilibrium. In a heavy ion collision, even if partons reach some initial local thermal 
equilibrium at early times, it is far from clear that the long wavelength modes of the 
chiral order parameter will be in equilibrium as chiral symmetry breaking occurs. 

In section 3.1, we will describe the phenomenological consequences of the existence 
of a large region of space in which the order parameter is displaced from the vacuum 
direction. Because of the explicit chiral symmetry breaking, oscillates about the 
cr direction. In later sections of this chapter, we will discuss a mechanism by which 
these long wavelength pion oscillations could arise.S To model the dynamics of the 
chiral order parameter in a far from thermal equilibrium phase transition, we will 
consider quenching in the 0(4) linear sigma model. We will argue, and present nu- 
merical evidence, that in the period immediately following a quench long wavelength 
modes of the pion field are amplified. This could have dramatic phenomenological 
consequences in heavy ion collisions. Simulations of classical evolution following a 
quench are a long way from real heavy ion collisions. A number of authors have 
discussed improvements and modifications — relaxing the quench approximation; in- 
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eluding effects of expansion; including quantum effects. Nevertheless, at this time 
a quantitative theoretical calculation is impossible. The biggest obstacle is that the 
appropriate initial conditions (initial meaning just before the transition occurs) on 
the long wavelength modes of the order parameter are not known. The physics of this 
chapter rests on much less solid theoretical foundations than that of the last chapter. 
Away from equilibrium, we cannot use renormalization group ideas or universality, 
and our conclusions are necessarily qualitative. If seen, long wavelength pion oscilla- 
tions and the consequent large fluctuations in the neutral to charged pion ratio would 
be a definitive signature of an out of equilibrium chiral transition. It remains an 
experimental question whether or not these phenomena occur in heavy ion collisions. 



3.1. Misalignment of the Chiral Condensate 

Among the most interesting speculations regarding ultra-high energy hadronic 
or heavy nucleus collisions is the idea that regions in which the chiral condensate 
is misaligned may occur.0iiii'iSfl'0 In such misaligned regions, the chiral conden- 
sate points in a different direction from that favoured in the ground state. If we 
parametrize the condensate using the variables of the sigma model, misaligned vac- 
uum regions are places where the four-component field = (a, vf), that in the ground 
state takes the value (w,0), is instead partially aligned in the vf directions. Because 
chiral symmetry is explicitly broken by small quark masses, the pion mass is not zero, 
and therefore if a misaligned vacuum region were produced, the field would oscillate 
about the a direction. Thus, such regions of misaligned condensate have a nonzero 
pion field oscillating in time. They would produce clusters of pions bunched in ra- 
pidity with highly non-Gaussian charge distributions. A misaligned vacuum region 
starting with the field in the tti — plane would emit only charged pions (equally 
positive and negative, since the fields are real), while a misaligned vacuum region 
starting with the field pointing in the tts direction would emit only neutral pions. 

More generally if we define 

/ ^ "'° , (66) 

then each misaligned region will yield a cluster of pions with some fixed /. Different 
regions, of course, will have different values of /. We can calculate the probability 
distribution of the ratio / if we assume that all initial values on the 3-sphere are 
equally likely. We will make this assumption in order to get a simple analytical result 
although it may not be strictly true because the quark mass term selects a preferred a 
direction even at high temperature. Where starts will determine in which direction 
it ends up oscillating about the sigma direction. We define angles on the 3-sphere 
according to 

(cr, TTs, TTi, 712) = (cos ^, sin COS sin ^ sin (yj COST/, sin ^ sin (yj sin?/) . (67) 
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Then, the ratio / is given by / = cos^ (p. Under the assumption that all initial values 
on the 3-sphere are equally probable, the probability distribution 'P{f) is determined 

by 



/ p[f)df = — df] dOsin^e ^ d^sin^ (68) 

J/l Tl^ Jo Jo Jarccos{y^) 

and turns out to be simply^ 

V{f) = \r''' . (69) 

Equivalently, the probability that / < /i is given by vTi- If large regions in which (p 
is oscillating coherently develop in a relativistic heavy ion collision, then there should 
be clusters of pions with small relative momentum in which / is constant, and the 
values / takes in different such regions should be distributed according to (|69D. 

As one application of (^), we note that the probability that the neutral pion frac- 
tion / is less than .01 is 0.1! This is a graphic illustration of how different (|69D is from 
what one would expect if individual pions were independently randomly distributed 
in isospin space. It also makes Centauro events in which less than 1% of the outgoing 
particles from a collision are neutral pions seem much less surprising than they first 
appeared. The analysis of the Centauro data is difficult for several reasons. Most 
important of these is the limitation imposed by small statistics. Also, if a Centauro 
event occurs too high above the detector, so many secondary photons will be pro- 
duced that the event will not be recognized as a Centauro. Third, in a Centauro event 
all of the particles from the collision strike a small region of the detector and it is 
impossible to isolate the central rapidity region. The balloon borne JACEE detectors 
do not have the second and third problems, but they are only exposed during brief 
balloon flights and therefore have very small statistics. When relativistic heavy ion 
collisions in a laboratory colliding beam facility occur at high enough energies that 
there is a low baryon number density central rapidity region, all of the difficulties of 
the cosmic ray experiments will rapidly be overcome. That will be the time to look for 
correlated clusters of pions, and to look for a distribution like (|69l). In the following 
sections, we propose and explore a concrete mechanism by which such phenomena 
might arise in heavy ion collisions in which the chiral order parameter is far from 
thermal equilibrium. We then discuss signatures in heavy ion collision experiments. 



3.2. Emergence of Long Wavelength Pion Oscillations After a Quench 

In Chapter 2, we considered the equilibrium phase structure of QCD. We argued 
that QCD with two massless quark flavours probably undergoes a second-order chiral 
phase transition. For many purposes it is a good approximation to treat the u and 
d quarks as approximately massless, and we expect that real QCD has a smooth but 

■^This result appeaxa^|t£j Jiajze^ been first known in 1981,@ and has been rediscovered a number of 
times since then.§@§SIS 
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perhaps rapid transition as a function of decreasing temperature from small intrinsic 
to large spontaneous chiral symmetry breaking. At first sight it might appear that 
a second-order phase transition is especially favourable for the development of large 
regions of misaligned vacuum. Indeed the long-lived, long-wavelength critical fluctua- 
tions which provide the classic signature of a second-order transition are such regions. 
Unfortunately we saw that light quark masses, even though they are formally much 
smaller than intrinsic QCD scales, spoil this possibility. The pion masses, or more 
precisely the inverse correlation length in the pion channel, are almost certainly not 
small compared to the transition temperature. As a result the misaligned regions are 
modest affairs at best even near the critical temperature. They almost certainly do 
not contain sufficient energy to radiate large numbers of pions. 

In this Chapter, we will consider an idealization that is in some ways opposite to 
that of thermal equilibrium, that is the occurence of a sudden quench from high to 
low temperatures, in which the (a, vf) fields are suddenly removed from contact with 
a high temperature heat bath and subsequently evolve mechanically. We shall show 
that long wavelength fluctuations of the light fields (the pions in QCD, which would 
be massless if not for the quark masses) can develop following a quench from some 
temperature T > Tc to T = 0. The long wavelength modes of the pion fields are 
unstable and grow relative to the short wavelength modes. 

In a real heavy ion collision, the phase transition must proceed by some process 
in between an equilibrium phase transition in which the temperature decreases arbi- 
trarily slowly and a quench in which thermal fluctuation ceases instantaneously. Our 
goal is to understand the dynamics of the long wavelength modes of the chiral order 
parameter in a quench in the hope that the qualitative behaviour in this model is 
representative of processes occurring in real heavy ion collisions in which the long 
wavelength modes are far from equilibrium, although not quenched. 

We shall use the linear sigma model of Gell-Mann and Levyil to describe the low 
energy interactions of pions at T = after a quench: 

C = j d^x\^ d'4>^d,K - \{<rK - v' f + ffcr I , (70) 

where is a four-component vector in internal space with entries (a, tt). Here, A, f, 
and H oc are to be thought of as parameters in the low energy effective theory 
obtained after integrating out heavy degrees of freedom. Unlike in the equilibrium case 
where universality allowed us to make quantitative predictions, at zero temperature 
we must make a non-universal choice of model Lagrangian to obtain equations of 
motion for the chiral order parameter. For now, we shall treat ([70|) as it stands as a 
classical field theory, since the phenomenon we are attempting to address is basically 
classical. We shall be dealing with energies of order 200 MeV or less, so that neglect 
of heavier fields seems reasonable. Since at zero temperature the a is heavy — it is 
a feature at about 600 MeV which is perhaps too broad to be called a resonance — 
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one might consider integrating it out. This would yield the nonlinear sigma model. 
However, we must implement initial conditions appropriate to a disordered state in 
which the tt and a fields are equivalent up to small effects due to explicit chiral 
symmetry breaking. Therefore, we must keep the sigma and use ([701). 

To model a quench, we begin at a temperature well above T^. The typical config- 
urations have short correlation lengths and (0) ~ 0. ((0) ^ because H ^ 0.) One 
then takes the temperature instantaneously to zero. The equilibrium configuration 
is an ordered state with the field aligned in the a direction throughout space, but 
this is not the configuration in which the system finds itself. The actual, disordered 
configuration then evolves according to the Lorentz invariant zero temperature equa- 
tions of motion obtained by varying ([70|) . Quenching in magnet models has been 
much studied in condensed matter physics.0 However in that context it is usually 
appropriate to use diffusive equations of motion, because the magnet is always in sig- 
nificant contact with other light modes {e.g., phonons). For this reason the condensed 
matter literature we are aware of does not directly apply to our problem. 

We will discuss numerical simulations of quenching to zero temperature in the 
linear sigma model with an explicit symmetry breaking term Ha which makes the 
pions massive. Turok and SpergelH have considered this scenario with no explicit 
symmetry breaking term as a cosmological model for large scale structure formation in 
the early universe. They find a scaling solution in which the size of correlated domains 
grows without bound at the speed of light. This is easily understood qualitatively, 
and is peculiar to the case where = 0. The field is at different points on the 
vacuum manifold in regions of the universe which have not been in causal contact 
after the quench. Because the vacuum manifold is degenerate in potential energy, 
the system evolves to reduce gradient energy. As time continues, larger and larger 
regions come into causal contact, align, and the size of correlated domains grow at 
the speed of light. When the 0{N) symmetry is explicitly broken, however, we do 
not expect a scaling solution. The Ha term tilts the potential, the vacuum manifold 
is not degenerate, and in a time of order the scalar field in all regions (whether 
in causal contact or not) will be oscillating about the sigma direction, and the physics 
cannot be described in terms of differently aligned domains coming into alignment. 

3.2.1. A Numerical Simulation 

We now describe numerical simulations of quenching in the linear sigma model 
which were reported in Ref . ||^ . We choose and randomly independently on each 
site of a cubic lattice. This means that the lattice spacing a represents the correlation 
length in the initial configuration. In the initial conditions, is disordered, the 0{N) 
symmetry is not spontaneously broken, and the lattice spacing a represents the vr and 
a correlation lengths which are approximately degenerate. In the simulation whose 
results are shown in Fig. 11, we chose and randomly from Gaussian distributions 
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centred around (p = ^ = and with (0^)^^^ = v/2 and (0^)^''^ = v. The three 
parameters v, H, and A in ([70D determine m^, mo-, and f^, = (0|cr|0) according to 



A(0|a|0)((0|a|0)" - v') - H = , (71) 
(0|(t|0) 

Note that f.,^ = {0\a\0) > v for H ^ 0. In interpreting our results we must remember 
that while in the code we are free to choose the energy scale by setting the lattice 
spacing a = 1, a actually represents the initial correlation length. In choosing the 
parameters for the simulation shown in Fig. 11, we assumed that a = (200 MeV)~^ 
and then chose v = 87.4 MeV = 0.4372 a-\ H = (119 MeV)^ = 0.2107 a'^, and 
A = 20.0 so that = 92.5 MeV, = 135 MeV, and = 600 MeV. 

With parameters and initial conditions chosen, we evolved the initial configuration 
according to the equations of motion using a standard finite difference, staggered 
leapfrog scheme. The simulation was performed on a 64^ lattice and the time step 
was a/10. After every two time steps, we computed the spatial fourier transform of 
0. For each component of 0, we computed the power 



(73) 

and averaged over all modes with k = \k\ in bins of width 0.057a ^. = (64a)'^ 
is the lattice volume. Fig. 11 shows the power in modes of one component of the 
pion field with spatial wave vectors of several different magnitudes as a function 
of time. The curves all start at approximately the same value at t = because the 
initial power spectrum is white since we chose independently at each lattice site. 
For a lattice with an infinite number of sites, with the normalization used in ([73|) 
the initial value of the power spectrum is related to the mean square of the initial 
distribution by P{k,0) = (1/A^)(0^), with = 4. This means that all the curves 
in Fig. 11 begin at t = at about 0.01 in lattice units. The behaviour of the low 
momentum pion modes is striking. While the initial power spectrum is white and, 
as ergodicity arguments would predict, the system at late times is approaching an 
equilibrium configuration in which the equipartition theorem holds, at intermediate 
times of order several times the low momentum pion modes are oscillating with 
large amphtudes. 

We verified that finite size effects are not important even for the longest wavelength 
{ka = 0.20) mode shown in Fig. 11 by checking that the behaviour of the ka = 0.31 
mode is the same in a 32^ and a 64'^ box. The longer the wavelength of a mode, 
the more that mode is amplified. This holds up to the longest wavelength mode for 
which we can trust the results of the simulation, namely the ka = 0.20 mode which 
corresponds to a wavelength A = 27r/A; of about half the lattice size. 



44 




Fig. 11. The time evolution of the power spectrum of the pion field at several spatial wavelengths 
after a quench at t = 0. The parameters of the potential and the initial configurations for and (f) 
are described in the text. Time is measured in lattice units, and since a = (200 MeV)^"'^ the unit 
of time is approximately 1 fm. The figure shows the time evolution of the angular averaged power 
spectrum (defined in the text) of one of the three components of the pion field. The curves plotted 
are (top to bottom in the figure) the average power in the modes in the momentum bins centred at 
ka = 0.20, 0.26, 0.31, 0.37, 0.48, 0.60, 0.71, 0.94, 1.16, 1.39, and 1.84. The initial power spectrum 
is white and all the curves start sd t — at approximately 0.01 in lattice units. The power in the 
longest wavelength pion modes is amplified by a factor of order 1000 relative to that in the shortest 
wavelength modes, which is not amplified at all. From Ref. [Q; extended to later times here than 
in Ref. g. 



45 



The behaviour of only one component of the pion field is shown in Fig. 11. The 
other two look qualitatively the same. It should be noted that the exact height of 
the peaks in the curves depend on the specific initial conditions. If the simulation is 
run with the same initial distributions for and but with a different seed for the 
random number generator, the heights of the peaks change, and the relative sizes of 
the peaks in the three different pion directions change. The qualitative features of 
Fig. 7 — the growth of long wavelength modes of the pion field — do not depend on 
the specific realization of the initial conditions. Amplification occurs for a range of 
choices of the parameters in the potential including physically motivated values. We 
will discuss the effects of varying the choice of initial distributions below. 

These numerical simulations have been repeated by Gavin, Gocksch and Pisarski,0 
who report that they have replicated the power spectra of Fig. 11. There has been 
some confusion in the literature, however. A number of authors have tried to describe 
the result seen in Fig. 11 in terms of a correlation length. In thermal equilibrium, this 
would be reasonable, as the power spectra would be completely specified by a single 
length — the correlation length. The essence of the phenomena of Fig. 11 is, however, 
that the spectra are nonthermal. They cannot be described by a single length scale, 
and computing a correlation length is meaningless. The classical simulation of Fig. 
11 ceases to be physically relevant at late times, when instead of classical fields one 
has noninteracting pions. Fig. 11 has been extended formally to late times, however, 
to show how at late times in the classical simulation the energy is equipartitioned 
and thermal equilibrium is reached. Hence, comparing intermediate and late times in 
the figure is instructive. At intermediate times, short wavelength modes have about 
the same amplitude as they do at late times in thermal equilibrium, while long wave- 
length modes have been amplified. It is impossible to describe this behaviour in terms 
of a single length scale. For a visual image, think of long wavelength ocean swells 
with short wavelength chop in addition; perhaps it is easier just to look at Fig. 11. 
The essential qualitative phenomenon whose cause and consequences we now discuss 
is that long wavelength pion oscillations have been excited — the short wavelength 
modes have not been damped out. This picture is more complicated than the idealized 
"smooth" region of disoriented chiral condensate described in the previous section. 

Although no meaningful correlation length or domain size can be defined, it is 
possible to describe the physics using the time dependent correlation functions 

Cix,t) = Jd'y^ix + y,my) (74) 

for each component of 0. For each component of 0, the power P{k,t) is the spatial 
Fourier transform of this correlation function. Since the essence of the physics we 
are describing is the amplification of long wavelength power, a description using the 
power spectrum as in Fig. 11 is the natural one, but the same information is of 
course encoded in the correlation function. Another reason why a description of the 
physics using P{k, t) is more natural than one using C{x, t) is that it is closer to what 



46 



experimentalists measure — they measure the momenta of particles not positions in 
space. Long wavelength oscillations seen as amplification of P{k, t) at some small 
magnitude k become pions with similar momenta (in the k direction) which are close 
to each other in an experimentalist's detector. 

3.2.2. Understanding the Results 

The emergence of long wavelength pion oscillations in the numerical simulations is 
a striking qualitative phenomenon. We now discuss a simple qualitative explanation 
for the result. Whenever one has spontaneous breaking of a continuous global sym- 
metry, massless Nambu-Goldstone bosons occur. The masslessness of these modes 
occurs because of a cancellation: 

= - fx^ + \<f , (75) 

where the second term arises from interaction with a condensate whose expectation 
value (0)^ satisfies (0)^ = /i^/A = in the ground state. However following 
a quench the condensate starts with its average at rather than v, and generally 
oscillates before settling to its final value. Whenever (0)^ < v^, will be negative, 
and sufficiently long wavelengths fluctuations in the Nambu-Goldstone boson field will 
grow exponentially. The same mechanism will work, though less efficiently, if one has 
only an approximate symmetry and approximate Nambu-Goldstone modes, or for 
that matter other modes whose effective is 'accidentally' pumped negative by 
their interactions with the condensate. 

Let us compare the growth of long wavelength modes found numerically with 
expectations based on the mechanism just discussed. Suppose the potential V in ([70|) 
were simply V{(f)) = (m^/2)(/)"0Q. Then, the equations of motion would be linear, and 
modes with different spatial wave vector k would be uncoupled. Each curve in Fig. 

11 would be a sinusoid with period vr/ \/ + P and constant amplitude. (The power 
spectrum, being quadratic in the fields, oscillates with one half the period of the 

fields.) The period of the oscillations in Fig. 11 is indeed given by vr/ \J + P, but 
the amplitudes are far from constant. This behaviour can be understood qualitatively 
by approximating 0^ in the nonlinear term in the equation of motion by its spatial 
average: 

0>„(f,t)~(0,0")(t) . (76) 

Note that the Hartree or mean field approximation (^) is not controlled. We make it 
only to obtain some understanding of the qualitative behaviour of the system. Using 
([76|) and doing the spatial fourier transform, the equation of motion for the pion field 
becomes ^ 

^7r{k,t) = -mlfJk,t)n{k,t) (77) 
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Fig. 12. Time evolution ol the spatially averaged {(f)'^) for the same simulation whose results are 
shown in Figure 1. The horizontal line is at (0^) — . From Ref. [Q. 

where the time dependent "mass" is given by 

ml^^{k,t) = -\v' ^k^ ^\{<\>^){t) (78) 

and where k = \k\. 

Fig. 12 shows the time evolution of (0^) in the same simulation whose results are 
shown in Fig. 11. In the initial conditions, (p is Gaussian distributed with {(f)'^) < v'^. 
Therefore, for a range of wavevectors with k less than some critical value, fnljj < 
and the long wavelength modes of the pion field start growing exponentially. (0^) 
grows and then executes damped oscillations about its ground state value (0|(7|0)^. 
A wave vector k mode of the pion field is unstable and grows exponentially whenever 
(0^) < v"^ — k'^/X. Modes with /c^ > Xv"^ can never be unstable. The k = mode 
is unstable during the periods of time when the (0^) curve in Fig. 12 is below 
f^. Since (0^) is oscillating about (0|o"|0)^ > f^, after some time the oscillations 
have damped enough that (0^) never drops below f^, and from that time on all 
modes are always stable and oscillatory. In general, longer wavelength modes are 
unstable for more and for longer intervals of time than shorter wavelength modes as 
(0^) oscillates. Also, m^jj- is more negative and more growth occurs for modes with 
smaller k. Making the approximation ( [76|) and thus using ( [77|) cannot be expected to 
completely reproduce the effects of the nonlinear term in (^) which is local in position 
space. Nevertheless, it predicts that the long wavelength modes of the pion field go 
through periods of exponential growth and predicts that the longer the wavelength the 
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greater the amplification, and therefore gives us some understanding of the behaviour 
of these modes in Fig. 11. In an equihbrium phase transition, exphcit symmetry 
breaking keeps the correlation length at Tc finite and, in the case of QCD, too short 
to be of interest. After a quench, on the other hand, arbitrarily long wavelength 
modes of the pion field are amplified even though the pion mass is nonzero. 

What about the sigma modes? In Ref. [|] the angular averaged power spectrum 
for the sigma field is plotted, as was done in Fig. 11 for one component of the pion 
field. For cr^ < 17^/3, the effective mass squared for the k = mode is negative. This 
condition is far less likely to be satisfied than {(j)"^) < v^. Hence, while some growth 
of the low momentum a modes is possible, they do not grow nearly as much as the 
low momentum pion modes do.@ 

At late times (0^) oscillates with small enough amplitude that it never goes below 
f ^ and consequently no modes are ever unstable. Hence, if we make the approxi- 
mation ( [76D we would expect that at late times each mode in Fig. 11 would con- 
tinue to oscillate with approximately constant amplitude, with the longer wavelength 
modes maintaining the large amplitudes acquired during their exponential growth 
spurts. That this is not what is seen in Fig. 11 reflects the effects we neglected 
in making the approximation (|76D. Because the modes are in fact coupled and the 
equations of motion are actually nonlinear, ergodicity arguments suggest that even- 
tually equipartition should apply. If the energy is equally divided among modes, 
then (amplitude)^ (m^ + k'^) should be approximately constant. Hence, the power in 
a mode at late times should be approximately constant in time and inversely pro- 
portional to {m\ + /c^), and this is indeed the behaviour at late times in Fig. 11. 
It is reasonable that longer wavelength modes take longer to decrease in amplitude 
than shorter wavelength modes. As the system heads toward equipartition, power 
has to flow from low momentum modes to higher momentum modes. This plausi- 
bly happens more quickly at higher momentum because higher momentum modes 
are coupled to a larger number of modes, or, in the continuum, a larger volume of 
momentum space. Hence, we end up with the striking phenomena of Fig. 11: while 
the system begins and ends with small n{k,t), in between the long wavelength pion 
modes are dramatically amplified. 

The late time behaviour of the simulation in Fig. 11 is not physical. When the 
expansion and cooling of the plasma is taken into account, there will be a time after 
which the energy density is low enough that the description in terms of classical fields 
no longer makes sense. After this time, one has individual pions flying off towards 
the detector. From Fig. 11, we see that if this occurs sometime in a broad range of 
times of about several tens of fm, it is plausible that at this time long wavelength 
modes of the pion field are greatly amplified compared to modes with higher \k\. 

It is important to note that we have only used the Hartree approximation (|76D as 
a qualitative guide. The emergence of long wavelength pion oscillations apparent in 
Fig. 11 occurs in a simulation of the strongly coupled nonlinear equations of motion. 
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We have already seen that the Hartree approximation leads to incorrect conclusions 
about the late time behaviour in Fig. 11. In fact, the Hartree approximation can 
lead one astray in another way if it is taken too literally. A number of authors have 
notedSS that (0^) of Fig. 12 is less than v"^ only at early times. Indeed, in the 
Hartree approximation, after a time of order ~ 1/3 fm, (0^) is oscillating about 
its vacuum value.Sil This suggests that the growth of long wavelength modes should 
cease after a time of this order. However, from Fig. 11 we see that such growth 
occurs for much longer times — until t ~ 30 fm. Thus, in this respect, the Hartree 
approximation fails to give a qualitative explanation of the phenomena observed in 
the full nonlinear equations. The possibility that the Hartree approximation can fail 
in this way was noted by Boyanovsky et al.^ Fig. 11 is evidence that it does fail in 
this way. 

In section 3.1, we discussed the phenomenological consequences of long wavelength 
oscillations of the pion field. Regions in which the oscillations are in the tts direc- 
tion become clusters of neutral pions, while regions in which the oscillations are in the 
TTi —TT2 plane become clusters of charged pions. Implicit in this description was the as- 
sumption that the latter clusters include pions with both positive and negative electric 
charge, and are approximately neutral. We can now justify this assumption.^ While 
the total electric charge must be conserved, there is no conservation law forbidding 
the emergence of regions of net positive and negative charge. We must determine 
whether long wavelength oscillations of the electric charge density p'^ occur. The 
electromagnetic current due to the charged pions is given by 

= - 7r2(9/,7ri (79) 

so the charge density is 

= 7ri7r2 - 712^1 . (80) 

It is easy to verify that the equations of motion obtained by varying (^) do indeed 
conserve j^. Because is conserved, we expect it to be more difficult to excite long 
wavelength charge oscillations than to excite long wavelength oscillations of tti, tc2, 
or 7r3. If we write (vri,7r2) as a complex field rexp(i6') then p^ = We have 

seen that there are periods of time when the pion field is exponentially unstable. 
This instability which manifests itself dramatically in Fig. 11 results in growth of 
the magnitude of tt. Thus, it leads to growth of r but does not tend to change 9. 
Equivalently, the instability tends to kick vf radially away from vf = 0, and results 
in oscillations with vf oscillating approximately through the origin with ^ ~ rather 
than with tt "orbiting" the origin with 6' 7^ 0. We have verified that the dramatic 
phenomena of Fig. 11 are not manifest in the electric charge density pe by computing 
the spatial fourier transform of in the simulation of Fig. 11. First, we have verified 
that the total charge in the 64'^ box is indeed conserved. Second, we found that the 
low k modes of p'^ do not grow dramatically. Long wavelength oscillations of tti, H2, 
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and TTs are dramatically amplified after a quench, while long wavelength oscillations 
of the conserved electric charge density pe are relatively unaffected. 



3.2.3. Other Scenarios 

Bjorken, Kowalski, and Taylor have proposed that regions of disoriented chiral 
condensate may arise in proton-proton collisions. ilSil They suggest that after a col- 
lision, a "cold" region throughout which the condensate is misaligned in a similar 
direction in isospin space is surrounded by a "hot" region of outgoing particles — 
hence the name Baked Alaska for the scenario. The MiniMax experimentli22l at Fer- 
milab (J. D. Bjorken and C. C. Taylor, co-spokespersons) is currently taking data to 
search for such phenomena. The scenario we have described is somewhat different, 
in that while long wavelength oscillations of the chiral condensate arise, we do not 
envision "smooth" regions of disoriented condensate. The long wavelength sloshing, 
although dramatic because it is completely different than what is obtained in thermal 
equilibrium, nevertheless has short wavelength oscillations superimposed on top of it. 
The emergence of long wavelength pion oscillations as described in this article requires 
that a large volume of space (certainly several fm across) is first in a high temper- 
ature phase in which the condensate is disordered {n.h disordered not disoriented). 
The volume has to be large enough that it makes sense to discuss long wavelength 
modes of the chiral order parameter. In this volume, the plasma must then experience 
a transition to the low temperature phase which occurs in such a fashion that the 
long wavelength modes of the order parameter do not stay in thermal equilibrium. 
We suggest that this physics may be relevant in the central rapidity region of a suffi- 
ciently energetic heavy ion collision, but it is not clear that it is possible to use this 
sort of language to describe a proton-proton collision. Although the Baked Alaska 
scenario is somewhat different than the one we are discussing, it has similarities. It is 
certainly far from thermal equilibrium. Furthermore, BjorkenE2ll has suggested that 
the processes at the boundary between the "hot" and "cold" regions of the Baked 
Alaska can be described as a quench. Finally, and most important, although the set- 
tings are different — hadron-hadron collisions vs. heavy ion collisions — at the end 
of the day long wavelength pion oscillations lead to large fluctuations in the neutral 
to charged pion ratio in either case. 

The idea that coherent pion emission from a classical pion field occurs in collision 
experiments has a long history,li°li3ll2lilS and solutions to the classical equations 



of motion have been studied in other contexts also.E^ In addition to the ideas we 
have described to this point, another possibility which has been discussed in the 
literature is that when there are a lot of pions in a small volume the fact that they 



are bosons can play a role — bose condensation can occur,E2a and the wave function 
must be symmetrized.0 This physics can have the effect of broadening the probability 
distribution for / from a narrow Gaussian, although it does not lead to a l/v7 
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distribution. Pehaps this could be associated with Centauro events. Greiner, Gong 
and MiilleJiS^ have argued that for these effects to be significant, pion densities larger 
than those likely to occur at RHIC are necessary, and that in addition if such large 
densities were to occur, a description in terms of a bose gas of pions would not be 
appropriate. There have also been recent proposals that classical pion fields could 



arise in hadron-hadron collisions " "^ H ^'^^ and in the leading (high-rapidity) regions of 
heavy ion collisions 
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It may turn out that the amplification of long wavelength modes described in 
this article may not be the appropriate explanation for the Centauro events, which 
may turn out to be hadron-hadron or hadron-nucleus collisions in which a large vol- 
ume in which the chiral order parameter in its high temperature phase is never 
formed. Related ideas like Bjorken, Kowalski, and Taylor's Baked Alaska scenario 
for the formation of a region of disoriented chiral condensate, or some of the other 
ideas just mentioned may turn out to be more appropriate descriptions of Centauro 
events, and the interesting phenomena which we hope Bjorken, Taylor, et al. discover 
experimentally^ in proton-proton collisions. 

In studying the chiral phase transition, perhaps the role of the Centauro events 
will prove to be merely inspirational, and not of direct relevance. Whatever the 
inspiration, it is important to consider the qualitative phenomena which are possible 
when the plasma in the central rapidity region of a heavy ion collision goes from 
the high temperature phase in which the chiral order parameter is disordered to 
the low temperature phase in such a way that long wavelength modes of the order 
parameter do not stay in equilibrium. We have used the quench as an idealization 
which may capture the relevant non-equilibrium physics. We have seen numerical 
evidence (supported by qualitative understanding) that long wavelength oscillations of 
the order parameter in the pion directions may be greatly amplified, whereas this does 
not happen if thermal equilibrium is maintained. In the next sections, we describe 
various improvements to the treatment, and discuss experimental consequences. 



3.3. Initial Conditions, Expansion, and Quantum Effects 

The results shown in Fig. 11 do depend on the initial distributions for (f) and (f). 
If the initial conditions are such that the energy per unit volume in the simulation is 
much more than |f no striking growth of long wavelength modes occurs. The system 
quickly achieves a configuration in which the energy is equipartitioned among modes 
corresponding to thermal equilibrium at a temperature above T^. Thus, in order for 
long wavelength modes to grow as in Fig. 11, the energy of the fluctuations in in the 
initial conditions must be less than it would be in thermal equilibrium. Another way of 
saying this is that whereas in the initial condition for Fig. 11 (0) is in a configuration 
appropriate to a high temperature, (0^) and (0^) have magnitudes characteristic of a 
lower temperature. These admittedly ad hoc quench initial conditions were chosen in 
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the hope of capturing quahtative behaviour typical of more general nonequilibrium 
physics. We now discuss efforts to improve upon this treatment. 



Gavin and Miilleilii3 note that if upon starting with quench initial conditions, the 
condensate were to evolve from (0) ~ to (0) ~ v much more slowly than it does 
in the classical linear sigma model considered above, then the amplification of long 
wavelength pion oscillations would be much greater than in Fig. 11. They call this 
scenario annealing, although this is a somewhat different usage than in condensed 
matter physics. Annealing usually means slow cooling in which both condensate and 
fluctuations stay close to their thermal equilibrium configurations as the temperature 
is slowly reduced. In this standard annealing, since approximate thermal equilibrium 
is maintained, no growth of long wavelength modes would occur. Gavin and Miiller's 
scenario, in which the condensate evolves slowly (as in standard annealing) but in 
which nonequilibrium conditions which allow amplification of long wavelength modes 
are maintained nevertheless, is in a sense more out of equilibrium than a quench. It 
would be marvelous if this turned out to be the appropriate idealization for treating 
heavy ion collisions, but that may be too much to hope for. 

A much more ambitious project than the classical simulations described so far 
is to include the effects of quantum fluctuations (in addition to classical thermal 
fluctuations) on the dynamics of the order parameter. A number of groups have 
made progress in this direction. SSIlll In order to treat the quantum fluctuations, 
however, all groups have had to make linearizing approximations either by working 
in the Nambu-Jona-Lasinio model, ii^ or by using the Hartree approximationS or the 
large-N approximation.Ell Such approximations cannot be avoided at present, but we 
have already seen that they are fraught with risk. Using the Hartree approximation 
to study the classical simulations of Fig. 11 leads to the prediction that growth of 
long wavelength modes occurs for a much shorter time than is in fact the case. Thus, 
it is perhaps not surprising that in treatments of a quench in which quantum effects 



have been included, short correlation lengths have been found. tlil'u Also, as we have 
noted before, the essence of the phenomenon discovered in the classical simulations 
cannot be described simply as the development of a long correlation length. Long 
wavelength modes are amplified, while short wavelength oscillations remain. Further 
work is needed in order to understand quantum effects on the dynamics. 

A less ambitious improvement on the simulations of the previous section is to 
include the effects of expansion. Following Bjorken,0 we wish to implement the 
boost invariant longitudinal expansion characteristic of the central rapidity region of 
a sufficiently energetic heavy ion collision. The natural coordinates to use are the 
rapidity y of equation (^) and the proper time r of equation (p2D, instead of the 
time t and the longitudinal position z. This can be implemented by taking the linear 
sigma model equations of motion and performing the replacement 
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The first treatments to include expansion were 1 + 1 dimensional — they assumed 
the configuration was uniform in the transverse directions.ii'EllBI "^ l ^^^ fellljll^ By 
making the further simplifying assumption that is independent of rapidity, Blaizot 
and Krzywickiii'0 reduce the problem to a + 1 dimensional nonlinear dynamical 
system. Following the evolution of these equations is instructive for analyzing the 
physics of a large smooth disoriented chiral condensate as envisioned in the Baked 
Alaska scenario of Bjorken et ai, and is a good toy model in which calculations can 
be done without resorting to numerical work."' KrzywickiEli has considered how the 
dynamics of this simple system can be affected by coupling to a bath representing all 
the other degrees of freedom in the system. Huang and Wang|iil have done a 1 + 1 
dimensional simulation in which and are random functions of rapidity with the 
same quench initial distributions that were used in the simulation of Fig. 11. They 
find that oscillations in the vr fields as a function of rapidity develop over distances 
in rapidity of up to 2 or 3. 

One reason why it is important to include expansion is that one can then ask 
whether it is possible to begin with initial conditions in thermal equilibrium at a 
high temperature and have something analogous to a quench occur as the system 
cools by expanding. This question has been addressed in numerical simulations by 
Asakawa, Huang, and Wangjiil and by Cooper et a/.S We will describe the work of 
both groups below, but let us first consider the qualitative answer to the question 
just posed. Both groups find that starting with thermal equilibrium, one stays close 
enough to thermal equilibrium that no significant growth of long wavelength modes 
occurs. Both groups also find, however, that if the long wavelength modes of the 
condensate are not in thermal equilibrium at the start of their simulations, then long 
wavelength pion oscillations can, in fact, be amplified. It would have been wonderful 
if simply adding expansion to the linear sigma model had been enough to turn thermal 
equilibrium into a quench, but perhaps this was too much to ask. Real heavy ion 
collisions are much more complicated than the linear sigma model, however. Indeed, 
using the sigma model to describe the dynamics cannot be valid for T ^ Tc. The 
unanswered question can be phrased as follows: In a heavy ion collision, at the time 
when it becomes appropriate to use the linear sigma model to describe the dynamics, 
are the long wavelength modes of the order parameter in thermal equilibrium with 
the short wavelength modes and the other degrees of freedom, or are they out of 
equilibrium? Nobody knows the answer. If they are out of equilibrium, then perhaps 
the quench discussed above or the quench with expansion to be discussed below yield 

Some author have used the further simplification gained by working in the nonhnear sigma 



model in which there is no a degree of freedom. In this article, we are interested in the amplification 
of long wavelength pion oscillations from an initial state in which is disordered and the tt and a 
fields are equivalent up to small effects due to explicit chiral symmetry breaking. We have therefore 
had to use the linear sigma model. If, however, one is interested in studying the evolution of long 
wavelength pion oscillations at later times without regard to how they may have arisen at earlier 
times, then the nonlinear sigma model suffices. 
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a qualitative description of the physics. If that is the case, then long wavelength pion 
oscillations arise. 



Asakawa, Huang, and Wangti3 have simulated the classical equations of motion 
of the linear sigma model in 2 + 1 dimensions, where the two spatial dimensions 
are the transverse ones. They include the effect of longitudinal expansion by using 
(l/r)9/9r(r9/(9r) instead of d"^ /dt^ in the equations of motion, but they do not do 
a 3 + 1 dimensional simulation. They are therefore able to use a large 2-dimensional 
lattice, and can implement transverse expansion. They do so by choosing initial 
conditions which are vacuum outside a cylindrical region of transverse radius 5 fm, 
and with disordered with a correlation length of 0.5 fm and with (0) = inside 
the cylindrical region. They do two simulations. The first has (0^) large enough 
to approximate thermal initial conditions, and not much happens in this simulation. 
The second simulation, shown in Fig. 13, has (0^) smaller, as in a quench. The 
transverse expansion is apparent in the figure, as is the growth of long wavelength 
pion oscillations. Just as in Fig. 11, one obtains a configuration with prominent 
long wavelength oscillations (about 4 fm from peak to trough; hence about 8 fm in 
wavelength) superimposed with short wavelength oscillations which have not been 
damped out. As in the simulations of Ref. 0], the longest wavelengths which are 
amplified are comparable to the system size. 

The linear sigma model simulations of Cooper et a/.^ll are 3 + 1 dimensional, in- 
clude the effects of quantum fluctuations, include the time dependent renormalization 
of the couplings in the potential due to the changing quantum and thermal fluctua- 
tions, include boost invariant longitudinal expansion, and do not include transverse 
expansion. In order to include quantum effects, they work in a mean field approx- 
imation, and so risk missing some of the physics as we discussed in section 3.2.2. 
Because longitudinal expansion is included, the amplitudes of all the fields becomes 
small at very late times, the equations of motion linearize, and it becomes possible 
to define and measure particle numbers. Results from their simulations are shown in 
Fig. 14. They plot the mean number density n{p) and its fluctuations as a function 



of transverse momentum p at several proper times. (See Ref. ||112|| for definitions, 
but growth of n{p) for low p corresponds to growth of long wavelength oscillations.) 
The left column is for a simulation with thermal initial conditions with an initial tem- 
perature of T = 200 MeV > Tc. Because of the longitudinal expansion, the system 
cools through Tc, but long wavelength modes are not amplified. The right column 
is for a simulation with nonthermal initial conditions in which (a) has been given 
an initial nonzero time derivative. Cooper et al. note that using quench-like initial 
conditions (long wavelength modes in a thermal distribution at one temperature with 
short wavelength modes in a thermal distribution at a lower temperature) produces 
similar results to those shown in the right column. Long wavelength oscillations are 
dramatically amplified. It is impossible to define a domain size at late times in the 
right column, because multiple length scales have emerged. For nonthermal initial 
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Fig. 13. From Ref. |119|. Contour plot of a component of the tt field (7r2) at the initial proper time 
T = 1 fm and later at r = 5 fm in a 2 + 1 dimensional simulation of the linear sigma model including 
expansion with quench initial conditions done by Asakawa et al. The uniform grey around the 
outside of both figures represents TT2 = 0. In the lower panel, there are prominent features at longer 
wavelengths than in the initial conditions shown in the upper panel. For example, the structure 
centred at about (2, 0) is a region of positive tt2 peaking at about 85 MeV, and the structure centred 
at about (—1.5,1.5) is a region of negative 7r2 with a minimum at about — —65 MeV. These 
structures are not smooth, however. (The grey-scale in this black and white representation of the 



figure does not make all this apparent — the original in Ref. |ll£l is in colour, and the authors 
also have a three dimensional colour version which is the clearest way of presenting the results, 
but which is not clear in black and white.) Thus, prominent long wavelength oscillations with 
wavelengths of about 8 fm have been excited. Because of the cylindrical boundary conditions used 
to implement the transverse expansion, these results cannot simply be Fourier transformed, as was 
done in producing Fig. 11. However, the phenomena are clearly analogous. Long wavelength modes 
have been markedly amplified, but they cannot be described as smooth domains as they have short 
wavelength oscillations superimposed on them. 
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Fig. 14. From the final paper of Ref. |112|. Results from simulations by Cooper et al. of the linear 
sigma model in a mean field approximation including quantum and thermal fluctuations and boost 
invariant longitudinal expansion. For nonthermal initial conditions, long wavelength oscillations are 
amplified, while this does not occur for thermal initial conditions. See text for further discussion. 



conditions, amplification of long wavelength modes is manifest in the particle number 
distribution, while for thermal initial conditions, this does not happen. This result is 
in qualitative agreement with that of Asakawa et al.j^ and confirms the observation 
first made in Ref. 0]: if the long wavelength modes of the order parameter are not 
in thermal equilibrium as the chiral transition occurs, there is a simple mechanism 
which can lead to dramatic amplification of long wavelength pion oscillations. 

To this point, the only quantum effects we have discussed are the quantum cor- 
rections to the dynamics. There are several other inherently quantum mechanical 
questions, however. First, once a long wavelength oscillation has been created, at 
late times as the classical treatment is becoming inappropriate one might worry that 
there may be enough quantum scattering of pions {e.g n^ir^ tt+tt^) to wash out 
the effect of interest before the pions are far enough apart on their way to the de- 
tector that they no longer interact. This question is being studied by Huang and 
Suzuki,0 and their preliminary result is that such effects are small because pions, as 
Goldstone bosons, have only derivative couplings and so have small scattering rates 
at low momentum. 

Another very interesting question which has received much attention is: To what 
quantum state does a classical configuration in which the chiral condensate is mis- 
aligned and oscillating correspond? In other words, once the oscillating classical field 
becomes noninteracting quantum mechanical pions, what is the miantum state of 



those pions? Following Horn and Silver,!^ a number of authorstiJ'tilJtliJllilZI have consid 
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ered the unique quantum state with pions with zero total isospin. Computing the 
distribution of the number ratio / of neutral pions obtained in this state does indeed 
yield the (l/2)/~^/^ distribution. Greiner, Gong, and Miiller note that this state also 
has unusual correlations. Defining 

e., ^ - 1 (82) 

where a and P can be 0, +, or — , they find that for states with many pions, Coo = 
4/5, Co+ = Co- = -2/5 and C++ = C__ = C+_ = 1/5. Note that in ordinary 
intensity interference experiments for bosons coherent sources give C = and 
C 7^ is only possible for identical bosons. Here, in contrast, neutral and charged 
pions are each positively correlated, and anti-correlated with the other. It seems 
likely that these charge dependent correlations will be qualitatively similar for other 
quantum states which may correspond to classical regions in which the condensate is 
misaligned. While it is reasonable to consider states with small isospin, there is no 



reason to require 7 = 0. It is therefore reassuring that Cohen et a/.llH3 have shown that 
the (l/2)/~i/2 distribution is significantly affected only when / becomes larger than 
about 30% of the number of pions. Once one allows consideration of quantum states 
which do not have a definite particle number, there are many possibilities. Amado 



and Koganll2^llHJ have enumerated several, including the isospin projected coherent 
states of Botke, Scalapino, and SugaJisl and squeezed states.tnl Other authorsSliii 
have noted that mixed states described by density matrices are a possibility. The 
argument for the relevance of squeezed stateJi^ is an interesting one. In the Hartree 
approximation, each mode of the pion field satisfies the equations of motion ([77| ) for 
a harmonic oscillator with a time dependent mass. Replacing the classical modes by 
quantum harmonic oscillators, one obtains parametric quantum oscillators, and this 



is just the setting in which squeezed states naturally arise.li^ The argument does rely 
on using the Hartree approximation, which we already know to be problematic, and 
so it is not conclusive. The question of what quantum state corresponds to a region 
in which long wavelength amplifications of a classical field remains open. 

3.4- Signatures and Experiments 

The definitive signature of long wavelength pion oscillations is large fluctuations 
in /, the neutral pion fraction, and we discuss this below. However, a number of other 



signatures have also been proposed. Huang, Suzuki, and Wang!l^'tH3 have noted that 
in a region in which the chiral condensate is misaligned, the isovector part of the 
electromagnetic current is no longer dominated by the p° meson as in the normal 
vacuum — it couples both to p mesons and to their chiral partners the ai mesons. 
Huang et al. consider a large smooth region of disoriented vacuum, and average over 
all possible directions of disorientation, and find that the p peak at 770 MeV in the 
dilepton spectrum will be reduced by a factor of 2 and a broad bump in the dilepton 
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spectrum will emerge at the ai energy (1260 MeV). The uj meson, on the other hand, 
is a chiral singlet and so will not be affected. They therefore propose that seeing a 
reduction in the p peak in the dilepton spectrum while the uo peak remains unchanged 
would be a signal for the formation of disoriented chiral condensate. It is not clear, 
unfortunately, that this signal would be prominent in the more realistic situation in 
which there are long wavelength oscillations of the condensate, but not large smooth 
regions. 

A direct signal of the amplification of long wavelength pion oscillations would be 
simply to look for an enhancement of the number of pions produced at low momentum. 



without regard to their charge.tH3'E^ It is not clear that such a signal would be 



unambiguous in and of itself. However, if this signal were detected in coincidence 
with the detection of large fluctuations in /, then it might provide a estimate of 
the wavelengths of importance — excess production of pions below some transverse 
momentum pt would show that modes with wavelengths longer than 1/pt had been 
amplified. 



It is also worth looking for the correlations CaR defined in (|82|). In particular Jl^ the 



unusual charge dependence of the correlations (neutral and charged pions correlated 



with themselves, and anti-correlated with the other) may be detectable. GavinllHa has 
also considered distinguishing the effect of regions of disoriented chiral condensate 
on the correlation between like sign charged pions from that of standard interference 
effects, but concludes that this will be difficult. 

We now turn to the definitive signature of long wavelength oscillations of the chiral 
order parameter. If a particular k mode ends up oscillating in the ttq direction, then 
it will become neutral pions moving in the direction of k. Similarly, oscillations in 
the TTi — plane will become charged pions. Thus, in some regions of the detector, 
/ would be large, and in others it would be small. If this were the only physics 
occurring, / would be distributed like (l/2)/~^/^ as we have discussed, but of course 
there will be a large background due to all the other pions in the event. (These 
correspond very roughly to the short wavelength oscillations which, as we have seen, 
remain even as long wavelength oscillations are amplified. This correspondence is 
only rough, however, because at short wavelengths the linear sigma model cannot 
capture all of the physics.) 

What is needed is a detector which covers as much of 47r solid angle as possible, 
which counts both charged pions and photons from the decay of neutral pions, and 
which is segmented in pseudorapidity rj and azimuth Lp into fine enough segments that 
the number of photons and charged pions in each segment can be accurately measured. 
(Ideally, a theorist's detector would be segmented enough that the typical segment 
had either zero or one photon or charged pion in it.) One would use this detector 
to look at the central rapidity region of the collision, where almost all the charged 
particles are pions and where almost all the photons come from neutral pion decay, 
so that particle identification is not crucial. It would be very helpful to have enough 
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information about the transverse momenta px so that high momentum charged pions 
and photons can be cut from the data. (In looking for a signal, one of the things to 
play with would be where to put the pr cut, but somewhere around 200 MeV would 
seem reasonable.) The data, then, would consist of an event by event catalogue of 
the positions in (77, y?) of the photons and charged pions. There are many conceivable 
ways of analyzing the data to look for large fluctuations in /. Several are described 
below. 

Let us first make a very crude attempt to estimate the number of pions that might 
arise from the excitation of a long wavelength pion oscillation. Inspired by Fig. 13, 
consider a wave with wavelength A ~ 8 fm and with k = 27c/X ~ 150 MeV. The 
energy density of such an oscillation with amplitude A is 

e^{ml + k'^)A^ ~ 25 MeV/fm^ . (83) 



The number of pions, then, is 




, \' , -125 (,,^)(^) (84) 



where V is the volume of space in which this mode is oscillating with this amplitude. 
The uncertainties in this estimate are huge. There is no information in the simulations 
of Fig. 13 about the longitudinal extent of the region which is oscillating. However, 
if it extends in rapidity from -1 to 1, then at a time t = 5 fm this corresponds to an 
extent in 2; of 8 fm. So, in ( |5^ we have taken ~ (8 fm)A^. Note, however, that 
there is no reason why V cannot be larger than — is the volume over which the 
mode is excited, and is not in principle related to A. We do not really know what 
amplitudes can arise — A ~ 75 MeV in Fig. 13, and in general it seems reasonable 
that A ~ is possible. Gavin has also tried to make an estimate of the number of 
pions from an event like that in Fig. 13;llHi he comes up with 30-300 pions, which 
seems reasonable. Following Horn and Silver ,0 Gavin@ notes that the momentum 
distribution of pions emitted from a smooth region of disoriented condensate of radius 
R goes like exp(— p^i?^). Taking R = A/2, this suggests typical transverse momenta 
of about 50 MeV. 

A few words of warning about the estimates just given seem in order. A reader 
at this point may be crying out that we should do a better job of analyzing the 
simulations done by various authors, and do better simulations, to make the estimates 
of the typical transverse momenta and pion multiplicity more accurate. This would 
be misleading. The theoretical foundations of the work described in this chapter 
are very speculative. Using the linear sigma model to describe the dynamics seems 
reasonable for the long wavelength modes, but cannot be justified in a controlled way. 
Either neglecting quantum effects or including them in a Hartree approximation risks 
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missing some of the physics. Most important of all, we do not know what initial 
conditions to impose on the long wavelength modes of the chiral order parameter at 
the time when partonic language ceases to be appropriate and a description of the 
dynamics in terms of the linear sigma model becomes appropriate. Thus, estimates 
of the kind given above are necessarily very crude. Indeed, the best advice for an 
experimentalist is to ignore all attempts at quantitative estimates, and just try to 
look at what is out there. Furthermore, although there are fewer complications at 
RHIC than at current lower energy experiments, given the lack of theoretical certainty 
it would be crazy not to look for fluctuations in / in current experiments at CERN 
and at the AGS. This is being or will soon be done,^^! but no results have yet been 
reported. 

Let us now consider how future data might be analyzed. The first thing to do will 
be simply to plot the data on a lego plot {i.e. in the (r/, v^) plane) using one colour 
for charged particles and another for photons, and then simply look by eye to see 
what is there. One of the consequences of the amplification of long wavelength pion 
oscillations is that rare dramatic events in which regions of the lego plot are almost 
all neutral or almost all charged will be much less rare than in thermal equilibrium. 
As KrzywickiEli has emphasized, even if such dramatic events remain rare, {i.e. less 
rare than in equilibrium but rare nevertheless) they are worth waiting for and looking 
for. 

Another simple way of analyzing the data will be to count all the charged particles 
and all the photons in each event, compute / ~ ^^^^ — jk for each event, and plot a 

I 1 '^charged ~r^7/*^ i 1 

histogram of the /'s for different events.t23 As an example, the WA98 experimentli^ 
at CERN detects both photons and charged particles in a window of rapidity of width 
1.1 near central rapidity, over all 27i of azimuth. In a typical event, of order 600 — 800 
charged particles are accepted. The task will be to compare the / histogram from 
the data with the / histogram for "background" events — presumably Monte Carlo 
generated — which do not contain any unusual fluctuations in / and so have a Gaus- 
sian / distribution with width proportional to the mean value of (n^/2 + ncharged)~^''^- 
Long wavelength oscillations of the chiral order parameter would lead to a broadening 
of this histogram. Clearly, understanding the / distribution of the background events 
is crucial. There is reason for optimism, however. It seems likely that if a Monte 
Carlo event generator predicts, say, more charged pions than actually occur, then 
it will also predict more neutral pions than actually occur, while the / distribution 
will be unaffected. This suggests that / distributions for background events obtained 
from event generators will be much more robust and accurate than, say, distribu- 
tions for the total multiplicity. This must of course be studied and tested further. 
Preliminary work suggests that a population of events given by taking Monte Carlo 
generated background events, superimposing 150 pions with momenta distributed like 
exp(— (p/50 MeV)^) in the centre of mass frame, uniformly distributed in azimuth, 
and with / distributed like (l/2)/~^/^, and then imposing WA98 acceptances, would 
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yield an / histogram which is clearly detectably broader than that of a population 
of background events. It might also be helpful to look for broadening of the / his- 
togram as a control parameter is varied which changes the energy of the collision. 
The beam energy is the obvious choice, but if that is not possible, looking at events 
with varying impact parameter may prove useful. At RHIC, the analysis analogous 
to the one sketched above would be to average over all charged particles and photons 
in a rapidity bin of some width, say 1, rather than to average over all particles in an 
entire event. 

It is important to realize, however, that averaging over azimuth is pessimistic. 
It is entirely possible that in a given event (or in a given region of rapidity in an 
event) long wavelength modes corresponding to one direction in azimuth will end 
up oscillating in one direction in vf space, while modes corresponding to a different 
direction in azimuth oscillate in a different direction in vr space. Stated in another 
way, it is clear that the pions from the long wavelength 712 oscillation in Fig. 13 
will not be uniformly distributed in azimuth. Thus, the generic signature of long 
wavelength pion oscillations is not large fluctuations of / which are uniform in cp, but 
rather large fluctuations in / which depend both on ip and rj. Thus, the data should 
be binned in bins of varying extents Arj and Aip, and / histograms should be done 
for varying bin sizes. Again, for each bin size, the / histogram from the data must 
be compared with that from background events. 



Another way of organizing the data analysis goes as follows.li^ Bin the data in such 
small bins that there are, say, about 10 particles per bin. Then for a bin centred at 
(?7, (p) define /(?7, (p) = f — {f) . (In a real experiment, once experimental acceptances 
are taken into account, (/) will not be precisely 1/3.) Because / is defined using very 
small bins, a histogram of / would be broad due to fluctuations even in the absence 
of interesting effects, and so would not be useful. However, compute the correlation 
function 

C{Ar^, A^) = (/(r/i, ^.0/(^2, ¥^2)) (85) 

where the averaging (...) is done over all bins with \rii — ri2\ = Arj and — (y92| = Aip. 
Seeing structure in this correlation function would suggest the appropriate scales in 
rapidity and azimuth on which to bin the data to see unusual / histograms. One could 
also do this analysis averaging over an ensemble of events. If all the events in the 
ensemble are selected from a single statistical distribution which is the sum of some 
signal distribution and some noise distribution, then the signal seen in C [Arj, Aip) 
would grow like the square root of the number of events in the ensemble. In this case, 
computing C{Ar], Aip) using an ensemble of events would be very valuable. However, 
if events are not selected from the same distribution, then averaging over an ensemble 
could wash out interesting effects. In particular, using an ensemble of events could 
miss seeing structures in which the interesting scales in rapidity and azimuth are 
different in different events, and could miss detecting effects which only occur rarely. 
Thus, the correlation function C [Arj, Aip) should be computed both event by event 
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and for ensembles of events. 

There are many ways to look for the fluctuations in the number ratio of neutral to 
charged pions which characterize long wavelength oscillations of the chiral condensate. 
We look forward to the day when the choice among analyses becomes data driven. 



3.5. Concluding Remarks 

Several different groupfllliiliii have now done complementary simulations in which 
different physics has been included and different assumptions have been made. All 
have found that long wavelength pion oscillations can be dramatically amplified if 
the chiral transition is out of thermal equilibrium. This behaviour can be understood 
at least qualitatively in terms of a simple mechanism: as the system evolves from a 
disordered configuration toward a low temperature ordered configuration, nonequi- 
librium dynamics of the condensate can give long wavelength pion modes a negative 
effective mass squared, making them unstable to growth. The abrupt transition, or 
quench, considered in section 3.2.1 is an idealization. In a heavy ion collision, the 
transition presumably proceeds by a mechanism which is in some sense between the 
quench idealization and the thermal equilibrium idealization. There remains room 
for improved simulations, complementary to all that have been done so far. A 3 + 1 
dimensional classical simulation with boost invariant longitudinal expansion done on 
a large enough lattice that transverse expansion can be included as was done in 2 + 1 
dimensions by Asakawa et al. would be very interesting. Furthermore, if the lattice is 
large enough that the simulation can be run to a late enough time that the classical 
field amplitudes linearize, then a direct computation of the particle number would 
be possible. Such a simulation would complement that of Cooper et al. because 
although it would not include quantum fluctuations, it would not require a mean 
field approximation and would include transverse expansion. However, none of these 
simulations can come close to describing a real heavy ion collision. First of all, the 
linear sigma model cannot capture all the physics, although it is likely a reasonable 
approximation for the long wavelength modes of the chiral order parameter. Perhaps 
the largest uncertainty arises from the fact that we do not know what initial con- 
ditions to impose on the a and tt fields at the time at which a linear sigma model 
description of the dynamics becomes appropriate. This is crucial, because simulations 
have shown that expansion beginning with thermal initial conditions is not sufficient 
to obtain the growth of long wavelength modes. It is necessary that at the time when 
the linear sigma model description becomes appropriate, the long wavelength modes 
of the order parameter are not in equilibrium with the short wavelength modes. 

The conclusions to take away from this chapter are essentially qualitative. Con- 
sider a heavy ion collision which is energetic enough that there is a central rapidity 
region of high energy density and low baryon number. If, as this region cools through 
the chiral transition, long wavelength modes do not stay in equilibrium, then there is 
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a robust mechanism which can lead to amphfication of long wavelength pion oscilla- 
tions. This cannot happen in thermal equilibrium. Long wavelength pion oscillations 
would be observed by detecting regions in rapidity-azimuth space in which the frac- 
tion of neutral pions is unusual. The pions from different long wavelength oscillations 
will have neutral pion fractions / with an inverse square root distribution. This will 
of course be masked by all the other pions in the events, but the signal to noise ratio 
can be increased by focussing on low pt pions. In contrast to the equilibrium physics 
of chapter 2 of this article in which we were able to use the notion of universality, 
the nonequilibrium physics which may lead to long wavelength pion oscillations is 
speculative. Therefore, this subject will only become mature when it becomes exper- 
imentally driven. A dedicated experiment to look for related phenomena is underway 
at Fermilab.ll2^ Experimentalists are looking for the phenomena described in this ar- 
ticle in current generation heavy ion experiments and will look for them in colliding 
beam experiments at RHIC and the LHC, where sufficient energy should be available. 
Detection of unusual fluctuations of / in the central rapidity region of a heavy ion 
collision would be a dramatic and definitive signature of an out of equilibrium chiral 
transition. The ball is in the experimentalists' court, and we wish them well. 
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